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We apply the theory of high-order harmonic generation by low- frequency laser fields in the strong 
field approximation to the study of the spatial and temporal coherence properties of the harmonics. 
We discuss the role of dynamically induced phases of the atomic polarization in determining the 
optimal phase matching conditions and angular distributions of harmonics. We demonstrate that the 
phase matching and the spatial coherence can be controlled by changing the focusing parameters 
of the fundamental laser beam. Then we present a detailed study of the temporal and spectral 
properties of harmonics. We discuss how the focusing conditions influence the individual harmonic 
spectra and time profiles, and how the intensity dependence of the dynamically induced phase leads 
to a chirp of the harmonic frequency. This phase modulation can be used to control the temporal 
and spectral properties of the harmonic radiation. Temporally, the harmonic chirped pulse can 
be recompressed to very small durations. Spectrally, chirping of the fundamental beam may be 
employed to compensate for the dynamically induced chirp and to control the individual harmonic 
spectrum. Finally, we discuss the short pulse effects, in particular nonadiabatic phenomena and the 
possibility of generating attosecond pulses. 
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I. INTRODUCTION 



A. Short history of high harmonic generation 



During the recent years high-order harmonic generation (HG) has become one of the major topics of super intense 
laser-atom physics. Generally speaking, high harmonics are generated when a short, intense laser pulse interacts 
with matter. Although HG has been in the course of recent years mainly studied in atomic gases, it has also been 
investigated in ions, molecules, atom clusters and solids. Apart from its fascinating fundamental aspects, HG has 
become one of the most promising ways of producing short-pulse coherent radiation in the XUV range. HG has 
already been a subject of several review articles |p4,109,123|. One can point out the following milestones in the short 
history of this subject: 



• First observations High harmonic generation is an entirely nonlinear and non-perturbative process. The spec- 
trum of high harmonics is characterized by a fall-off for the few low order harmonics, followed by an exte nded 
plateau, and by a rapid cut-off. The first experimental observations of the plateau were accomplished by [106] 
and §| in the end of 1980's. 

• Plateau extention Most of the early work has concentrated on the extension of the plateau, i.e. generation of 
harmonics of higher frequency and shorter wavelength [105,^l|. By focusing short-pulse terawatt lasers in rare 



gas jets, wavelengths as short as 7.4 nm (143rd harmonic of a 1053 nm Nd-Glass laser [117]), 6.7 nm (37th 
harmonic of a 248 nm KrF laser [ 12C| ]) and 4.7 nm (169th harmonic of an 800 nm Ti-Sapphire laser j3^|) have 
been obtained. Very recently, with ultra short intense infrared pulses, it has be come p ossible to generate XUV 
radiation extending to the water window (below the carbon K-edge at 4.4nm) [144,33. 



1 A breakthrough in the theoretical understanding of HG process in low frequency laser fields 
I who have shown that the cut-off position in the harmonic spectrum follows the universal 



• Simple man's theort 
was initiated by |8C 

law I p + 3C/ p , where I p is the ionization potential, whereas U p — e 2 £ 2 / Amuj 2 , is the ponderomotive potential, 
i.e. the mean kinetic energy acquired by an electron oscillating in the laser field. Here e is the electron charge, 
m is its mass, and £ and to are the laser electric field and its frequency, respectively. Pretty soon an explanation 
of this universal fact in the framework of "simple man's theory" was found [^ , p8[ . According to this theory, 
harmonic generation occurs in the following manner: first the electron tunnels out from the nucleus through 
the Coulomb energy barrier modified by the presence of the (relatively slowly varying) electric field of the laser. 
It then undergoes oscillations in the field, during which the influence of the Coulomb force from the nucleus 
is practically negligible. Finally, if the electron comes back to the vicinity of the nucleus, it may recombine 
back to the ground state, thus producing a photon of energy I p plus the kinetic energy acquired during the 
oscillatory motion. According to classical mechanics, the maximal kinetic energy that the electron can gain is 
indeed ~ 3U p . A fully quantum mechanical theory, that is based on strong field approximation and that recovers 
the "simple man's theory" , was formulated soon after ^,^9| . 

• Ellipticity studies The "simple man's theory" leads to the immediate consequence that harmonic generation 
in elliptically polarized fields should be strongly suppressed, since the electron released from the nucleus in 
such fields practically never comes back, and thus cannot recombine |58|,[39| . Several groups have demonstrated 
this effect p9|j46] , p9[ | , and have since then performed systematic experimental |30,165,9 166 141] and theoretical 
[p2|,p|j2^] studies of the polarization properties of harmonics generated by elliptically polarized fields. 

• Optimization and control Progress in experimental techniques and theoretical understanding has stimulated 
numerous studies of optimization and control of HG depending on various parameters of the laser and the active 
medium. These studies involved among others: 
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Optimization of laser parameters These studies concer n fo r instance laser polarization (discussed above), 
pulse duration or wavelength dependence [jl5 75 3^Jl| , 132 ] . Although typically infrared (Nd-Glass o r Ti - 



Sapphire) lase rs ar e used, generation by shorter wavelength intense KrF lasers is also very efficient ( [ 120 1 , 
for theory see [136])- 



Generation by multicolored fields Harmonic generation in combined laser fields of t wo frequ e ncies was 
studied in the context o f (i) enhancement of conversion efficiency (for theory see jj4,122 150 77, ll(f| , for 
experiment cf. [ 161 , 1 12 1 ) , (ii) access to new frequencies and tunability, if one of the fields is tunable (for 
theory see 



2|, for experiments |53|]6l| ]), and (iii) the control of HG process in general J67[ . 

Optimization of the generating medium First of all, optimization with respect to atomic gases was studied 
Other active medi a apart f ro m no ble gases have been use d to generate harmonics: ions ( thes e involve 



104], atomic 



ionized noble gas atoms [ 138 158 , It . 12C ] , and alkaline ions [|2|, [159| ) , molecular gases 
clusters p9[ | etc. At this point it is worth adding that harmonic generation from solid targets and lase r 
induced plasmas has also been intensively studied in the recent years (for theory see [63,125,100 130 131], 
for experiments cf. ]3l|,f74], pM| , pTl[ ) . 

Optimization and characteristics of spatial and temporal properties Those studies concern spa tial, temporal, 
and spectral properties of harmonic radiation, and in particular their coherence properties [132]; they are 
closely related to the subject of this review and will be discussed separatel y bel ow. Other examples of such 
studies involve spa tial c ontrol of HG using s patia lly dependent ellipticity [107], control of ph ase m atching 
conditions for low 108 and high harmonics |l34| , role of ionisation and defocusing effects 109 47]. 

Other control schemes. Other control schemes of har monic ge neration have been proposed that involve for 
example the coherent superposition of atomic states [ 163, 137 1 



• Applications High harmonics provide a very promising source of coherent XUV radiation, with numerous ap- 
plications in various areas of physics. In particular, applications in atomic physics are reviewed by Jl6| , p8[ . 
Harmonics have already been used for solid state spectroscopy p5| , and plasma diagnostics 152 1 . Further 
applications that employ directly coherence properties of harmonics will be discussed in this review. 

• Attosecond physics Future applications of high harmonics will presumably involve attosecond physics, i.e. the 
physics of generation, control, detection and application of sub-femtosecond laser pulses. Two types of proposals 
how to reach the subfemtosecond limit have been p ut forw ard over the last few years: those that rel y on pha se 
locking between consecutive harmonics [F8]]^M|39|37,160|, and those that concern single harmonics [ 139, 132| . 



B. Spatio-temporal characteristics of high harmonics 



1. Experiments 

Both from the fundamental and practical points of view it is very important to know and understand the spatial 
and temporal coherence properties of high harmonics. Informations on the spatial coherence of the beam and on its 
focusability, its spectrum and time profile are of direct interest for applications. But they also help in understanding 
the physics of the process, since there are many possible causes for distortion of the spatial and temporal profiles, and 
their interpretation implies a rather refined and deep study of the problem. 

The spatial distribution of the har moni c emission has been investigated by several groups in various experimental 
conditions. Peatross and Meyerhofer |l!3j ] used a 1 /im 1 ps Nd-Glass laser loosely focused (f/70) into a very diluted 
gaseous media (1 Torr) in order to get rid of distortions induced by phase matching and propagation in the medium. 
The far-field distributions of the harmonics (11 to 41) generated in heavy rare gases were found to be quite distorted, 
with pedestals surrounding a narrow central peak. These wings were attributed to the rapid variation of the harmonic 
dipole phase w ith t he laser intensity. 



Tisch et al. [151] studied high-order harmonics (71 to 111) generated by a similar laser, focused (f/50) in 10 Torr 
of helium. Complex spatial distributions are found for harmonics in the plateau region of the spectrum. However, in 
the cutoff, the measured angular distributions narrow to approximately that predicted by lowest-order perturbation 
theory. The broad distributions with numerous substructures observed in the plateau are attributed to the influence 
of ionization, and in particular of the free electrons, on phase matching. 

The influence of ionization on spatial profiles has also been investigated exp erimentally by L'Huillier and Balcou 1 9^] 
for low-order harmonics in xenon, and, more recently, by Wahlstrom et al. [159] for harmonics generated by rare-gas- ike 



ions. Generally speaking, ionization induces a significant distortion of the harmonic profiles, thus complicating their 



interpretation. In a recent Letter [133], we presented results of an experimental study of spatial profiles of harmonics 
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generated by a 140 fs Cr:LiSrAlF6 (Cr:LiSAF) laser system. Thanks to this very short pulse duration, it was possible 
to expose the medium to high intensities while keeping a weak degree of ionization. Under certain conditions, the 
resulting harm onic profiles were found to be very smooth, Gaussian to near flat-top, without substructure. 



In Ref. [135 1, we present systematic experimental studies of harmonic angular distributions, investigating the influ- 
ence of different parameters, such as laser intensity, nonlinear order, nature of the gas and position of the laser focus 
relative to the generating medium. We show that, when the laser is focused before the atomic medium, harmonics 
with regular spatial profiles can be generated with reasonable conversion efficiency. Their divergence does not depend 
directly on the nonlinear order, the intensity or even the nature of the generating gas, but rather on the region of 
the spectrum the considered harmonic belongs to, which is determined by the combination of the three preceding 
elements. When the focus is drawn closer to the medium, the distributions get increasingly distorted, becoming 
annular with a significant divergence for a focus right within or after the jet. 

A first endeavour to measure the degree of spatial coherence of the harmonic radiation has been recently made by 
Ditmire et al. p^ j with a Young two-slit experiment. They investigated how the coherence between two points chosen 
to be located symmetrically relative to the propagation axis depends on the degree of ionization of the medium. 
Finally, very recently far-field interference pattern created by overlaping in space tw o beams of the 13th harmonic, 



generated independently at different places in a xenon gas jet was observed [160 167 1. 

The experimental studies of temporal and spectral properties of high-order harmonics have also been carried out in 
different experimental conditions. Temporal profiles of low harmonics generated by relativel y lon g pulses (several tens 
of ps) were measured using a VUV streak camera by Faldon et al. ]57j and Starc zewski et al. |145| ] . In order to measure 
the duration of harmonic pulses in the femtosecond regime, Schins et al. ]140[ developed a cross-correlation method 
in which they ionize helium atoms by combining two pulses: the fundamental (800 nm, 150 fs) from a Ti:Sapphire 
laser, and its 21-st harmonic (38 nm). These pulses generate characteristic electron spectra whose sidebands scale as 
the cross-correlation function which can be mapped out by varying the delay between the two pulses. Using variants 
of this method Bouhal et al. p7j , and Glover et al. |54l were able to measure the duration of 21-st to 27-th harmonics 
within a sub-picosecond accuracy. For instance, in Ref. |27j] for a fundamental pulse of 190 fs FWHM, durations of 
100±30 fs and 150 ± 30 were found for the FWHM of 21-st and 27-th harmonic, respectively. 

Conc ernin g the spectral properties of the individual harmonics, the blue shift due to ionization has been rep orted 



in Ref. [158], whereas spectral properties of harmonics generated by chirped pulses have been discussed in Ref. [ 168 



2. Theory 

The theoretical description of spatial distributions and temporal profiles of harmonics requires to combine a reliable 
single atom theory that describes the nonlinear atomic response to the fundamental field, with a propagation code 
that accounts for phase matching, dispersion etc. Peatross et al. have studied the spatial profiles of low-order 



harmonics in the loosely focused regime [114 115|. Muffet et al. [ 1 10 1 modeled the results of Ref. [151], and showed 



that, depending on the focusing conditions, substructures could be either due to ionization or to resonances in the 



intensity d epen dence of the atomic phase. Temporal profiles of low-order harmonics were discussed in Ref. |57,145| 



Rae et al. [128 1 performed calculations outside the slowly varying envelope approximation by solving simultaneously 
the equations for the atomic dynamics and propagation, using a one dimensional approximation. Temporal and 
spectral profiles were studied in the strongly ionizing regime. 

In the series of papers p7| , p9[ , we have developed a single atom theory which is a quantum-mechanical version of 
the two-step model of Kulandcr et al. B4] and Corkum |3q] . This theory has been combined with the theory of HG 



by macroscopic media p5| to describe experimental results in a realistic manner. In a recent Letter |134| we have 
stressed the role of the dynamically induced phase of the atomic polarization in phase matching and propagation 
processes. In particular we have demonstrated the possibility of controlling the spatial and temporal coherence of the 
harmonics by changing the focusing conditions of the fundamental. We have performed numerical simulations of the 
angular distributions. Th e sim ulated profiles reproduced remarkably well the experimental trends and are thus used 
to interpret them in Ref. |135| . The role of the intensity dependent phase of atomic dipoles was elaborated in more 
detail in Ref. |)0) (see also |70|) . In the Ref. (!(]] we present a short review of the various consequences of the intensity 
dependent phase. 

C. The scope of the present review 

As already stated above, the knowledge of the coherence properties of high harmonics is of major importance both 
for applications and from the fundamental point of view. Although experimental and theoretical work has been 
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already devoted to this subject, a systematic study of the spatial and temporal coherence properties of harmonics is 
still lacking. In particular, the implications of the existence of a phase of the harmonic dipole have not been fully 
explored. The aim of this review is to present a detailed theoretical study of this important subject. The theoretical 
approach used in this paper is very well established, and has been confronted numerous times with experimental 
results with great success. Some of the results presented in this review have been published before and compared 
with experiments. Many of the results, however, are either new, or have been only reported in the PhD thesis of P. 



Salicrcs [ 132 1 . Nevertheless, we feci that very soon these results will find their experimental confirmation. 

To some extent this review has a character of a case study, i.e. we discuss here in great detail quantitative coherence 
properties of a specific high harmonic in specific conditions etc. It is important, however, to keep in mind that the 
presented results are not only qualitatively, but also quantitatively valid in more general cases. We hope that these 
results will turn out uscfull for anybody interested in applications of harmonics in general, and their extraordinary 
coherence properties in particular. 

The plan of the review is the following. In Chapter 2, we give an overview of theoretical approaches and describe 
our theoretical method: the single atom theory based on strong field approximation, and the propagation equations. 
Since our theory has been discussed in detail in other publications, we limit ourselves to present the final expressions 
that we use for calculations of the physical quantities. In Chapter 3 we discu ss the phase matching problem stressing 



the role of the dynamically induced phase of the atomic polarization [134 90 1. The combined effects of this phase 
and the phase of the fundamental beam depend on the atomic jet position relative to the focus. We investigate the 
influence of the jet position on the conversion efficiency. We show how the phase matching effects modify the cut-off 
law. 

Chapter 4 starts with a short section devoted to the general definitions of the degree of coherence and characteristics 
of partially coherent beams. We then concentrate on the emission profiles and quality of the wavefronts of harmonics, 
both in the near-field and in the far-field zones. We also present calculations of the degree of spatial coherence of the 
harmonics. 

In Chapter 5 we turn to the discussion of the temporal and spectral coherence. We show how the intensity 
dependence of the phase of the atomic polarization leads to a temporal modulation of the harmonic phase and 
to a chirp of its frequency. In both Chapters 4 and 5, we relate our theoret ical findings to experimental results. 



In particular, we use the parameters corresponding to experiments of Rcf. [ 135 ], and discuss systematically the 



dependence of the coherence properties of harmonics on the focus position of the fundamental. This parameter, as 



shown in Ref. [134], allows us to control the degree of coherence; optimal coherence properties are obtained when the 
fundamental is focused sufficiently before the atomic jet. In the final part of Chapter 5, we discuss the possibility of 
optimizing and/or controlling the temporal and spectral properties of harmonics by making use of the dynamically 
induced chirp: temporal compression with a grating pair and spectral compression by using a chirped fundamental 



pulse. These ideas are confronted with the recent experiment of Zhou et al. [168], and to the theoretical proposal of 



Kulander [139| 



In Chapter 6 we discuss future applications of harmonics with the special emphasis on their coherence properties. 
In particular we focus on applications in interferometry and on the short pulse effects. We discuss the possibility of 
generating and applying attosecond pulses. Finally, we conclude in Chapter 7. 

II. THEORY OF HARMONIC GENERATION IN MACROSCOPIC MEDIA 

The theory of harmonic generation in macroscopic media must necessarily contain two components: (i) a single atom 
theory that describes the response of an atom to the driving fundamental laser field, and (ii) a theory of propagation 
of the generated harmonics in the medium. In this Chapter we outline various approaches to describe these two 
components of the theory. 

A. Single atom theories 

The single atom theory should describe the single atom response to a time-varying field of arbitrary intensity, 
polarization and phase. In other words, it should allow to calculate the induced atomic polarization, or dipole 
acceleration, which then can be inserted as a source in the propagation (Maxwell) equations. In principle it is 
sufficient to describe the atomic response in the framework of the single active electron (SAE) approximation (cf. j32[; 



for a discussion of two-electron effects see for instance 87,55,145 ]). Also for relatively long laser pulses (of duration 



down to ~ 50 — 100 fs for Nd-Glass or Ti-Sapphire lasers) one can use the adiabatic approximation, i.e. calculate the 
atomic response for the field of constant intensity, and only at the end integrate the results over the "slowly" varying 
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envelope of the laser pulse. The discussion of the validity of the adiabatic approximation is presented in more detail 
in Chapter 5. 

There are essentially four methods that have been used to solve the problem of the single atom response : 

• Numerical methods These methods allow to solve the time dependent Schodinger equation (TDSE) describing an 
atom in the laser field. Since (at least in th e adi abatic case) the field oscillates periodically, one of the possible 
approaches is to use the Floquet analysis ( [ 1 19 1 ; f or a recent review see (6lJ), but the direct integration of the 
TDSE is far more often used (for a review see [ 33 , 123 1 ) . In ID such integration can be performed using either 
the finite element (Crank-Nicholson), or split operator techniques; in the context of harmonic generation it has 
been first used by the Rochester group pH , but then employed by many others as a test method. 

In 3D the numerical method has been initiated by Kulander ||sl| , [S2|| , who used a 2D finite element ("grid") 
method. Pretty soon it was realized that basis expansion methods that employ the symmetry of the problem 
(i.e. the spherical symmetr y o f a bare atom, or the cylindrical symmetry of an atom in the linearly polarized 
field) work much better f45| , |S5|] . Modern codes use typically expansions in angular momentum basis, and solve 
the coupled set of equations for the radial wave function using finite grid methods (cf. J80|), Sturmian expansions 
[H,[ll) or B-spline expansions [§l],U| . Most of those cod es ar e quite powerful and allow to calculate the atomic 
response directly without adiabatic approximation (cf. [139]). Unfortunately, they are quite time consuming, 
and it is therefore very hard to combine the results obtained from the numerical solutions of TDSE with the 
propagation codes. The reason is that the single atom response in the physically interesting regime is typically 
a rapidly varying function of the laser intensity and other laser parameters. The propagation codes thus require 
very detailed data from single atom codes. This problem becomes even more serious in the absence of cylindrical 
symmerty; real 3D numerical codes (such as the ones describing generation by elliptically polarized fields) have 
been developed only recently (cf. j^Jll,124]), and obviously are even more time and memory consuming. 

Nevertheless, many seminal results concerning harmonic generation has been obtained using direct numerical 



methods: from the first observation of the /„ 



-SU p law 



to the recent proposal of attosecond pulse generation 
139 1 . Particularly interesting are the contributions of the O xford-Im perial College group (for a review see ll23| ) 
that concern among others HG by s hort wavelength lasers [ 120, 136 1, pulse sh ape and blue shifting effects jL62| , 
role of strong ionization in HG [127], tempor al aspect s of harmonic emission [128,164], and the generation from 
the coherent superpositi on of a tomi c states | 163j , 137 [. The TDSE method has also been applied to molecules 
aligned in the laser field [169,79 118|. 



• Classical phase space avaraging method A lot of useful information about high harmonic generation processes 
can be gained from a purely classical analysis of the electron driven by the laser field. In order to mimic quantum 
dynamics, classical Newton equations are solved here for an ensemble of trajectories generated from an initial 
electron distribution in the phase space. This distribution is supposed to mimic the true quantum initial state 
of the system, so that averages over this distribution are analogs of quant um averages. Such approach has been 
developed in the context of HG by Maquet and his collaborators Jl9,155| (sec also fisfl). 



• Strong field approximation As already mentioned, the seminal paper on the I p + 3U P law [ pCfl stimulated the 
formulation of the "simple man's theory" |84|]38] 1 . Originally, this theory has been formulated as a mixture of 
quantum and classical elements: first, the tunneling of the electron out from the nucleus was described using 
the standard ADK ||[|4],[78| theory of tunneling ionization. The subsequent oscillations of the electron in the 
laser field were described using classical mechanics. Finally, electron recombination back to the ground state 
was calculated using the classical cross-section for the collision and the quantum mechanical recombination 
probability. 

A fully quantum mechanical theory, that recovers the "simple man's theory" in the semi-classical limit was 
formulated soon after |97],|8!|. This theory is based on the strong fi eld a pproximation (SFA) to the TDSE. It 
is a generalization of the Keldysh-Faisal-Reiss approximation | f7l| , p6 ,129|, applied to the problem of harmonic 
generation. It was for the first time formulated in the context of harmonics by Ehlotzky |52j| ; it is also strongly 
related to the, so called, Becker model of an atom with a zero range pseudo-potential interacting with the laser 
field (see below). 

In our formulation, the theory is based on the following assumptions: (i) it neglects all bound states of the 
electron in an atom with exception of the ground state; (ii) all the states in the ionization continuum are taken 
into account, but in their dynamics only the part of the Hamiltonian responsible for the oscillations of the free 
electron in the laser field is kept. Technically, we disregard all off-diagonal continuum-continuum transitions that 
change electronic velocity. With this two assumptions, the TDSE becomes exactly soluble, and the resulting 
solutions are valid provided U p > I p . In most of applications we treat the electronic states in the continuum as 
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Volkov plane waves, which additionally limits the validity of our approximation to electronic states with high 
kinetic energy, and thus to the description of the generation of high harmonics (with photon energy > I p ). It is 
worth stressing, however, that a heuristic scheme of accounting for Coulomb potential effects in the continuum 
within the framework of SFA was proposed by |(|| . 

Originally our method has been formulated for linearly polarized laser fields in the adiabatic (slowly varying 
intensity) approximation We have since then generalized it to elliptically polarized fields ^ , two-color 

fields |62], and the fields with periodically time-dependent polarization fl^ ], where all of those results were 
obtained in the adiabatic approximation. Finally, the method has been generalized to the fields with arbitrary 
time-dependence without adiabatic approximation (see Chapter 5). We have also performed intensive studies 
of the semiclassical approximation applied to our method in order to understand the harmonic emission in 
terms of semiclassical electron trajectories and Feynman path integrals |9(],[l(| . These studi es we re essential for 



understanding the role of the intensity dependent phase of the nonlinear atomic polarization [134]. They allowed 
also to understand the mechanism of electron trajectory selection in propagation, responsible for the generation 
of attosecond pulse trains |7|,|l7j|. The advantage of our method, apart from its very transparent physical sense, is 
that it gives partially analytic results, allowing rapid calculation of the very precise data required for propagation 
codes. Last, but not least, our method combined with propagation codes gives results in very good agreement 
with experiment; it has become the standard theoretical method of analysis of experimental data in the Saclay 
and Lund groups; it is also used by other groups |77| , |50| ] . 

• Pseudo-potential model Many important results in the theory of harmonic generation have been obtained by 
Becker and his collaborators who have solved exactly (and to a great extent analytically) the zero range pseudo- 
potential model j2(J . In this model the electron is bound to the nucleus via the potential 

V{r) = — <5(r)|-r, (1) 
Km or 

where m is the electron mass. This potential supports a single bound state with the energy —I p = — k 2 /2m. 
This model, ori gina lly formulated in the case of a linearly polarized field, was also extended to one-color [ p2[ 



and two-color [ 103 fields with arbitrary polarization. It w as also used to study the polarization properties 
of harmonics generated by elliptically polarized fields |102[] . As our SFA theory |8^^|, Becker's model may 
rigorously account for the ground state depletion |]23|| . Structures in the harmonic spectra were associated in this 
model to the above-threshold ionization channel closings pl[ | , rather than with quantum interferences between 
the contributions of different electronic trajectories Nevertheless, Becker's model leads practically to the 
same final formulas for the induced atomic dipole moment as our SFA theory, and to very similar results (the 
small discrepancies are caused by additional approximations used for numerical elaboration of final expressions; 
for detailed comparison of the tw o mo dels see p4|). Becker's model has also been used by several groups to 



analyze experimental data [105 54 112] 



B. Single atom response in the strong field approximation 



In this section we present explicit formulas describing the response of a single atom to the laser pulse in the s trong 
field approximation. Since the details of our version of the SFA can be found in the series of Refs. [37,39 134 9Cj,8fi, 
we limit ourselves to present the relevant expressions and to discuss their physical sense. 

Within our approach we obtain an approximate solution of the time dependent Schrodinger equation that describes 
an atom in the strong electric field of a laser of frequency to in the single active electron approximation. The 
knowledge of the time dependent wave function ^(i)) allows us to calculate the time-dependent dipole moment 
x(t) = (ty (t)\x\ty (t)) in the form of a generalized Landau-Dyhne formula (we use atomic units) 

f°° ( 7T \ ^ 2 

*(*) = * / dr ( — d*{p s - A(t)) exp (-iS(p s , t, r)) 

Jo \v + it/2J 

x £(t - T)-d(p s - A(t - t)) exp(-Tt) + c.c, (2) 

where v is a positive regularization constant, A(t) denotes the vector potential of the electromagnetic field, £(t) = 
—dA(t)dt = (Ei x cos(wi), sin(wt), 0) is the electric field (polarized elliptically, in general), whereas S(p,t,r) is 
the quasi-classical action, describing the motion of an electron moving in the laser field with a constant canonical 
momentum p, 
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.S'(/7./.ri= / dt" A P ? +I p \, (>!) 



with I p denoting the ionization potential. In expression (g) we have already performed (using the saddle point method) 
the integral over all possible values of the momentum p with which the electron is born in the continuum. For this 
reason the integral in Eq. (|J) extends only over the possible return times of the electron, i.e. the times it spends in 
the continuum between the moments of tunneling from the ground state to the continuum and recombination back 
to the ground state. The saddle point value of the momentum (that is at the same time a stationary point of the 
quasi-classical action) is 

p s =p s (t,r) = f dt" A(t")r. (4) 



Note the characteristic prefactor (y + ir/2)~ 3 / 2 in (|J) coming from the effect of quantum diffusion. It cuts off very 
efficiently the contributions from large t's and allows us to extend the integration range from to infinity. 

The field-free dipole transition element from the ground state to the continuum state characterized by the momen- 
tum p can be approximated by |S9| , p5f 

2 7 / 2 a 5 / 4 p 

d{p) = i — — — -, (5) 

IT [P + Oi)^ 

with a — 2I p , for the case of hydrogen-like atoms and transitions from s-states. 

Finally, T is the ionization rate from the ground state. In the framework of our theory it can be represented as 
twice the real part of the time-averaged complex decay rate 

f°° f 7T \ - 

lit) = dr £*{t) ■ d*(p s - A(t)) exp(-iS(f s , t, r)) 

Jo \v + it/2) 

x£(t-r)-d(p s -A(t-r)). (6) 

Note that both expressions (0) and (0) have the characteristic form of semiclassical expressions that can be analyzed 
in the spirit of Feynman path integral: they contain (from right to left) transition elements from the ground state to 
the continuum at t — t, propagator in the continuum proportional to the exponential of i times the quasi-classical 
action, and the final transition elements from the continuum to the ground state. Applying the saddle point technique 
to calculate the integral over r (and t if one calculates the corresponding Fourier components or time averages) , one 
can transform both expressions into the sums of contributions corresponding to quasi-classical electron trajectories, 
characterized by the moment when the electron is born in the continuum t s — r s , its canonical momentum p s (t s , t s ) (see 
Eq. (|J)), and the moment when it recombines t s Note, however, that due to the fact that we deal here with the 

tunneling process (i.e. passing through the classically forbidden region), these trajectories will in general be complex. 
Typically, only the trajectories with the shortest return times Re(r) contribute significantly to the expression (||); 
there are two such relevant trajectories with return times shorter than one period, i.e. < Re(ri) < Re(r2) < 2tt/lo. 

Note also that the dipole moment (|J) can be written in the form 

x{t) = x q e- iquJt - Tt + c.c. (7) 

q odd 

where x q denote Fourier components. They can be calculated either directly from Eq. (^|) using a Fast Fourier 
Transform, or analytically as discussed in Ref. ||. 

It is important to remember that both expressions (||) and (||) result from the single active electron approximation. 
Before inserting these expressions into the propagation equations, one has to account for the contributions of all active 
electrons, and replace Eqs. (||) and (||) by the total dipole moment and the total ionization rate that are given by the 
sums of the corresponding (independent) contributions of all active electrons. In the case of Helium (two s electrons 
in the ground state), this amounts to multiplying both expressions by the factor two. In the case of other noble 
gases (six p electrons in the ground state, two in each of the m = —1,0, 1 states) the procedure is more complex. 
Both expressions should be replaced by two times the sums of contributions of the given magnetic quantum number 
to = —1,0, 1; each of those contributions should be calculated replacing Eq. (||) by an appropriate field-free dipole 
matrix element describing the transition from the I = l,m = — 1,0, 1 states to the continuum. 

Fortunately, the dependence of the dipole moment and ionization rate on the details of the ground state wave 
function is rather weak, and typically reduces to an overall prefactor (p9Lp|, that determines the strength of the 
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dipole, but not the form of its intensity dependence. For these reasons, in most of the calculations for noble gases 
other than Helium, we still use the s-wave function to describe the ground state but multiply the results by an 
effective number of active electrons, n e i, that is equal to that ranges between 2 < n e ; ~ 4 < 6 for other noble gas 
atoms. Total ionization rates of Helium and Neon calculated with n e i = 2, ~ 4, respectively, agree very well with the 
ADK ionization rates 



C. Propagation theory 



In order to calculate the macroscopic response of the system, one has to solve the Maxwell equations for the fun- 
damental and harmonic fields. This can be done using the slowly varying envelope and paraxial approximations. The 
fundamentals of such an approach have been formulated by . Several groups have used simila r approaches to study 
the effects of phase matching, and to perform direct comparison of the theory with experiments 1 110 , 114 , 1 15 126 128 ]. 
To our knowledge, the most systematic studies of these sort have been so far realized by the Saclay-Livermore-Lund 
collaboration. 

In a series of papers we have studied propagation and phase matching effects in the context of the following problems: 
(i) phase matching enhancement in nonpertubative regime fj3| ; (ii) phase matching effects in tight focusing conditions 
|B6|,n ; (iii) shift of the observed cut-off position | p7[ ; (iv) density dependence of the harmonic generation efficiency 
Pf; (v) harmonic generation by elliptically polarized fields Q; (vi) harmonic generation by two-colored fields [ |G2| ; 
(vii) coherence control of harmonics by adapt ing t he focusing conditions [134]; (viii) influence of the experimental 
parameters on the harmonic emission profiles 135 1; (ix) generation of attosecond pulse trains fit]; (x) generation of 
attosec ond puls es by laser fields with time-dependent polarization jn| , (xi) interference of two overlapping harmonic 
beams [160 167 1 , and others. We have also formulated generalized phase matching conditions that take into account 
the intensity dependent phase of the induced atomic dipoles fl7| l . 



D. Macroscopic response 

In this section we present the Maxwell equations for the fundamental and harmonic fields used in our above 
mentioned studies. Using the slowly varying envelope and paraxial approximations, the propagation equations can be 
reduced to the form (here we use SI units) 

Vi^ 1 (r,t) + 2^fc; aEl n (r '* ) +2k%6k 1 (r,t)E 1 (r,t)=0 (8) 
oz 

BE (r t) 

V 2 ± E q (r, t) + 2ik f gg2 ! + 2k° g [Ak° g (z) + Sk q (r, t)]E q (f, t) 

= -^4p q NL (r,t), (9) 

where Ei(f,t), and E q (f,t) denote the slowly varying (complex) envelopes of the fundamental and harmonic fields 
respectively, k q = qto/c, whereas the rest of the symbols is explained below. The slow time dependence in the 
above equations accounts for the temporal profile of the fundamental field that enters Eq. (||) through the boundary 
condition for E\. The solutions of the propagation equations for given t have therefore to be integrated over t. 

The terms containing Ak g (r,t) describe dispersion effects due to the linear polarisability of atoms, and can be 
in fact neglected in the regime of parameters considered (low density). The terms proportional to 5k q (r,t) — 
—e 2 J\f e (r,t)/2mqcLu, with e denoting the electron charge, m- its mass, and Af e (r,t) the electronic density, describe 
the corrections to the index of refraction due to ionization; here the ionic part of those corrections is neglected. The 
electronic density is equal to the number of ionized atoms, i.e. 



Af e (r,t)=Af a (z) 



l-exp(-J TdE^f^Ddt' 



(10) 



where N a (z) is the initial density of the atomic jet, and T(\Ei(r,t')\) is the total ionization rate, which takes into 
account the contributions of all active electrons calculated from Eq. (J6j) using T = 2Re[/ j(t)dt/T] with T = 2n/u 
for an instantaneous and local value of the electric field envelope Ei(r, t'). Note that since T depends functionally on 
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Ei(f,t'), En. (|§) is a nonlinear integro-differential equation; it has to be solved first, and its solution is used then to 
solve Eq. (||). 

Finally, the Fourier components of the atomic polarization are given by 

P q NL (r, t) = 2Af a (z)x q (r, t)e^ (r ^ exp (- j T{\Ei{?, t')\) dt^j , (11) 

where x q (r,t) denote the harmonic components of the total atomic dipole moment, that includes the contibutions 
of all active electrons calculated from Eq. (|J) for a field (I-Eizl cos(u;i), \E\ y \ sin(wi),0). The factor of 2 arises from 
different conventions used in the definitions of P q L and x q . Finally, 4>i(r,t) represents the phase of the laser field 
envelope E± (r , i) , obtained by solving the propagation equation for the fundamental. 



III. PHASE MATCHING 



In this chapter, we study how the harmonic field builds up in the medium, and is influenced by the dynamically 
induced phase of the atom ic polari zation. In all the calculations presented in this work, we try to mimic the experi- 
mental conditions of Refs. [135,134]. The 825 nm wavelength laser is assumed to be Gaussian in space and time, with 
a 150 fs full width at half maximum (FWHM). The laser confocal parameter b is equal to 5 mm and the focus position 
is located at z = Q. The generating gas is neon and the atomic density profile is a Lorentzian function centered at z 
with a 0.8 mm FWHM, truncated at z ±0.8 mm. 



A. Source of the harmonic emission 



The generic intensity dependence of the dipole strength (i.e. absolute value squared) and phase for the 45th harmonic 
generated by a single neon atom is shown in Fig. |lj. At low intensity, when the harmonic is in the cutoff region, 
the dipole strength (dashed line) increases rather steeply with laser intensity and the phase (solid line) decreases 
linearly. When the harmonic enters the plateau region, the strength saturates and exhibits many interferences, while 
the phase decreases twice as fast as in the cutoff region, predominantly linearly but with superimposed oscillations. 
The regular behaviour of the dipole in the cutoff region is due to the existence of only one main electron trajectory 
leading to the emission of the considered harmonic. On the contrary, in the plateau, there are at least two relevant 
trajectories whose contributions interfere, resulting in a perturbed dipole. In this region the dominant contribution 
is typically from the trajectory with the longer return time ti (see section 2.2) - it is the action along this trajectory 
which determines the mean slope of the intensity dependence of the phase. 

One should mention, however, that the mean slope of the intensity dependence is not always a good measure of 
the phase variations in the macroscopic medium. The reason is that propagation and phase matching may lead to a 
single trajectory selection that depends on the relative position of the atomic jet with respect to the focus 0. For 
the case when the jet is before, or close to the focus, phase matching usually selects the dominant trajectory with 
the return time T2, and the intensity dependence of the phase after the selection follows the mean slope of the total 
phase. When the jet is after the focus, however, and when intensity in the jet is high enough, the trajectory with the 
return time T\ is selected, and the slope of the intensity dependence is much weaker than the mean slope of the total 
phase. The latter case will not be considered in the following; we shall discuss only the situations when the mean 
slope of the phase can be used to characterize the intensity dependence of the phase in the macroscopic medium. 

This intensity dependence of the dipole is quite universal for sufficiently high order harmonics, and depends weakly 
on harmonic number, laser wavelength, generating gas. These parameters only determine the position of the plateau- 
cutoff transition. In the following, we will thus concentrate on the study of the 45th harmonic generated in neon, 
keeping in mind that its behaviour is very general. 



B. Dynamically induced phase of the atomic polarization 

Harmonic generation is optimized when phase-matching is achieved, i.e. when the difference of phase between the 
generated field and the driving polarization A<I> — $ g — $ po / is minimized over the medium length, allowing an efficient 
energy transfer. In the case of plane waves and low order harmonics, this results in the well known phase matching 
condition on the wave vectors: Ak = k q — qk\ « 0. However, for the very high orders discussed here, the variation of 
the phase of the polarization is much more important than that of the induced field, so that matching phases amounts 
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to minimizing the phase variation of the driving polarization. One of the main causes of variation of the polarization 
in the medium is due to the rapid variation of the dipole phase with intensity. It induces a spatially and temporally 
dependent phase term in the medium, which influences the generation of the macroscopic field. 




FIG. 1. 45th single neon atom harmonic intensity (dashed line) and phase (solid line) as a function of the laser intensity. 

Let us consider this variation on the propagation axis at the maximum of the pulse temporal envelope for the 45th 
harmonic. It is induced by the variation of the intensity I(z) = Iq/(1 + Az 2 /b 2 ), and is shown by the short-dashed 
line in Fig. || for a peak intensity, i.e. at best focus and at the maximum of the pulse envelope, of 6 xlO 14 W/cm 2 . 
The variation of the dipole phase is less rapid outside the interval [—3 mm, +3 mm ], when the intensity on axis 
corresponds to the cutoff region (I(z) < 2.4 x 10 14 W/cm 2 , see Fig. [l]). The other important contribution to the 
polarization phase is a propagation term induced by the phase shift of the Gaussian fundamental field in the focus 
region, equal to — q tan~ 1 (2z/6), q denoting the process order. This function is shown in long-dashed line in Fig. 
^. There are other possible causes of variation of the polarization phase, such as atomic or electronic dispersion, but 
they are negligible in the conditions considered here. The total phase of the nonlinear polarization is represented in 
Fig. H by the solid line. In the region z < 0, the variations of both phases add, leading to a rapid decrease of the total 
phase. In the region z > 0, they have opposite signs, and almost compensate when the intensity on axis corresponds 
to the cutoff region. Consequently, phase matching strongly depends on the position of the medium relative to the 
laser focus. The best phase matching conditions on axis are those for which the phase variation of the polarization 
over the medium length (~ 1 mm) is minimal, i.e. when the laser is focused approximately 3 mm before the generating 
medium. Note that, at the minimum of the curve close to the focus, the superimposed oscillations are detrimental to 
a good phase matching. 
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FIG. 2. Phase of the polarization on the propagation axis (solid line). The long-dashed line indicates the term due to the 
propagation of the fundamental, and the short-dashed line the dipole phase for a peak intensity of 6 xlO 14 W/cm 2 . The laser 
propagates from the left to the right. 
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FIG. 3. Phase of the polarization for different radial positions relative to the propagation axis, for a peak intensity of 6 
xlO 14 W/cm 2 : r — (short-dashed), r — 5/im (solid), r = 10/im (long-dashed) and r — 15/im (dot-dashed). The dotted 
lines indicate the edges of a gas jet placed in z — —1mm, and the horizontal solid line, a trajectory r(z) that keeps the phase 
constant. 



So far, we have only considered phase matching on axis, which corresponds to centered harmonic profiles. However, 
good phase matching off axis can be realized in certain conditions. This is illustrated in Fig. ^ with the longitudinal 
variation of the polarization phase for different radial positions, from r =0 to 15 /an (relative to the propagation axis) 
for a peak intensity of 6 xlO W/cm 2 . Here, for clarity, the curves have been smoothed so that the superimposed 
oscillations in the plateau region do not appear. For a gas jet centered in z — — 1 mm, it is possible to minimize the 
longitudinal phase variation by moving from one curve to the other, i.e. by going off the propagation axis. Along these 
favored directions, the quick variation of the laser intensity on axis (cause of the rapid decrease of the polarization 
phase) is avoided by going off axis. In these conditions, the harmonic field can build up efficiently in the plateau 
region. Note that a method allowing the systematic study of these generalized phase-matching conditions has been 
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proposed in Ref. fL7[ . 
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FIG. 4. Phase of the polarization on the propagation axis. From the top to the bottom, the laser intensity I = 2, 3, 4, 5, 6 
xlO 14 W/cm 2 . 

If we now consider another peak intensity, the shape of the total phase is modified. This is illustrated on axis in 
Fig. |||for several peak intensities, from 2 to 6 xlO 14 W/cm 2 . As the intensity increases, the induced phase becomes 
more and more important in determining the total phase variation near the focus, which departs more and more from 
the arctangent term. 

The optimal phase matching position on axis is observed at different z depending on the peak intensity, since it 
always corresponds to the plateau-cutoff transition of the dipole (2.4 xlO 14 W/cm 2 ). Thus, for a given geometry, 
there will not be a static phase matching during the laser temporal envelope, but a continuous distortion of the build 
up pattern in the medium. This dynamic phase matching complicates the interpretation of the processes. 

C. Influence of the jet position on the conversion efficiency 




Jet/focus position (mm) 



FIG. 5. Conversion efficiency for the 45th harmonic as a function of the position of the center of the jet relative to the focus, 
for peak intensities ranging from 3 to 6 xlO W/cm . The laser temporal envelope is square with a 150 fs width. 
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FIG. 6. Conversion efficiency for the 45th harmonic as a function of the position of the center of the jet relative to the focus, 
for peak intensities ranging from 3 to 6 xlO 14 W/cm 2 . The dots show the influence of ionization at 6 xlO 14 W/cm 2 . The laser 
temporal envelope is Gaussian with a 150 fs FWHM. 



Using the numerical methods described in Chapter 2, we perform the propagation of the generated harmonic fields 
in the medium, considering for the moment square laser temporal envelopes, i.e. static phase matching. In Fig. |^, 
we study the variation of the conversion efficiency for the 45th harmonic generation as a function of the position z of 
the center of the atomic medium (relative to the laser focus placed in z — 0), for peak intensities ranging from 3 to 6 
xlO 14 W/cm 2 . The peak atomic density is 15 Torr. 

At low intensity, the curve presents only one maximum located in z = 1 mm. Increasing the intensity, this 
maximum splits into two lobes that become more and more separated. This is the confirmation of the two optimal 
phase matching positions described above. The positions of the maxima z > correspond precisely in Fig. M to the 
best phase matching positions on axis, i.e. to the plateau-cutoff transition of the dipole (2.4 xlO 14 W/cnr). The 
maxima occuring for negative z correspond to optimized phase matching positions off axis, as shown in Fig. |3J. Note 
that all the curves arc asymmetric compared to the focus position (z = 0), and that the conversion efficiency in z = 1 
mm is larger at the lowest intensity. However, this effect disappears when we consider harmonic generation by a 
Gaussian laser pulse, as shown in Fig. ^ for the same peak intensities as before. During the major part of the laser 
pulse, the intensity is in the cutoff region where the polarization amplitude drops, resulting in a lower conversion 
efficiency than for a square pulse. Except for this modification, the general behavior is the same as for square pulses, 
with enlarged peaks but similar positions and efficiencies. 

Note that if we take into account the ionization of the medium at 6 xlO 14 W/cm 2 (dots in Fig. ^|), we find a 
marginal influence on the number of photons except very close to the focus position. In the following, we will thus 
neglect ionization in a first step, and include it afterwards when its effects are not negligible. Wahlstrom et al. [158] 
have measured about 10 5 generated photons for the 45th harmonic in neon at 6 xlO 14 W/cm 2 . The difference of 
one order of magnitude with the results of our simulations can be explained, at least partly, by the more optimized 
conditions used in their experiment (higher pressure and longer confocal parameter). 
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FIG. 7. Intensity dependence of the conversion efficiency for the 45th harmonic for a jet position in z = (solid line) and 
z — 1 mm (dot-dashed line). Ionization is taken into account. The short-dashed line indicates the strength of the dipole 
moment (arb. units). 

D. Modified cut-off law 

The influence of the position of the focus relative to the gas jet on the harmonic conversion efficiency has an 
interesting consequence on the dependence of the harmonic yield as a function of intensity, for the different geometries. 
Fig. |?] presents the intensity dependence of the conversion efficiency for the 45th harmonic for a focusing at the center 
of the jet (solid line, z—0). The comparison with the strength of the dipole moment (short-dashed line) indicates two 
main consequences of the propagation: the rapid variations in the plateau region are smoothed out, and the change 
of slope indicating the plateau-cutoff transition is shifted to a higher intensity. This shift is the same for higher order 
harmonics, and thus implies that propagation decreases the extent of the plateau of the harmonic spectrum compared 
to the single atom response, from a photon energy of I p + 3.2U P to about I p + 2U p . This was observed experimentally 
by L'Huillier et al. p7| , and explained in terms of the variation of phase matching with intensity. 

If the laser was not focused in the jet, one would expect that the lower intensity experienced by the nonlinear 
medium would result in an even larger shift of the plateau-cutoff transition (recall that we are plotting the curves as 
a function of the peak intensity, i.e. at the focus and at the maximum of the pulse envelope). In Fig. [7j the intensity 
dependence for the case when the jet is located in z — 1 mm is shown in dot-dashed line. Amazingly, the shift of the 
plateau-cutoff transition to higher intensities is less important than for z = 0, corresponding to a cutoff law of about 
I p + 2.3U p . This is a direct consequence of the optimization of phase matching at low intensity for this particular 
position, as shown in Fig. || This effect is independent of the nonlinear order considered, and would indicate that 
the maximal extent of the plateau is not obtained for a focusing right into the jet, but rather slightly before it. 

IV. SPATIAL COHERENCE 
A. Definition 

We recall here some notions of the theory of partial coherence (see, for example, Born and Wolf pq|). The 
coherence of a beam is related to the correlation of the temporal fluctuations of the electromagnetic fields inside 
this beam. It is thus characterized by its mutual coherence function, defined for any two points inside the beam by: 
FiaCr) = (E\(t + r)-E|(t)), where E\, E2 are the complex amplitudes of the electric field in these two points, and the 
angular brackets denote an appropriate time average (here, over the harmonic pulse). The normalized form of the 
mutual coherence function is the complex degree of coherence: 
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(T) = r ia (r) = {E« + r)m)) (12) 

Vrn(o)r 22 (o) V(l^il 2 )(l^l 2 ) 

whose modulus is known as the degree of coherence. The temporal coherence is described by 711 (r), while the spatial 
coherence is described by 712(0). Note that the latter, despite its name, is related to the correlation in time of the 
fields emitted in two points. These quantities determine the ability of the fields to interfere, and can be measured 
in interferometry experiments: In Young's two-slit experiment, with both slits uniformly illuminated, |^12 (0) | at the 
slit positions is simply given by the fringe visibility, defined as V = (l ma x ~ I m in)/{I m ax + Imin) (where I max and 
Imin are the maximum and minimum intensities of the fringe pattern). The spatial coherence length of a beam at a 
given distance from its focus is defined as the length over which the degree of spatial coherence is larger than some 
prescribed value (between 0.5 and 0.9, depending on the authors, and on the coherence of their own source...). 

Another important aspect related to the coherence of a beam is the quality of its wavefront, an aspect that is often 
confused with the preceding description. A beam is said to be "diffraction limited" if the product of its spot size (at 
focus) and of its far-field spread (divergence) is of the order of the wavelength. This is realized when both the focal 
spot presents a reasonably regular amplitude variation, and the phase front at focus is very well behaved (typically 
plane). In particular, any distortion of the phase front will result in a larger (e.g. N times) angular spread, and the 



beam will be called "N times diffraction limited" 142 j 



In the following, we shall concentrate on the two focusing positions corresponding to well defined phase matching 
conditions at 6 xlO 14 W/cm 2 , namely z — 3 mm (on axis) and z = — 1 mm (off axis). For these extremal positions, 
on either side of the conversion efficiency curve (see Fig. |^), phase matching is mostly efficient close to the maximum 
of the laser temporal envelope, thus simplifying the study. Note that the main dependences of th e ha rmonic emission 



profiles (laser intensity, nonlinear order, jet/focus position) have been intensively studied in Ref. 135], and compared 
successfully with experimental data. We here focus on the coherence properties. 



B. Study of the spatial coherence: atomic jet after the focus 




FIG. 8. Normalized spatial profiles for the 45th harmonic at the maximum of the pulse for a jet position in z = 3mm, for 
intensities from 4 to 6 x 10 14 W/cm 2 . In solid lines are shown the profiles at the exit of the medium as a function of the radial 
coordinate, and in dashed lines, the far-field profiles (divergence). 

First we study the characteristics of the harmonic beam in the near-field, i.e. at the exit of the medium in z = 3.8 
mm (the half width of the jet is 0.8 mm) and at the maximum of the laser pulse. Fig. || presents in solid lines 
the harmonic profiles corresponding to different intensities, from 4 to 6 xlO 14 W/cm 2 (square pulses). The first two 
profiles are Gaussian, with 12 and 14 /xm radius in 1/e 2 respectively, while the third is super Gaussian with a 20 
/Mii radius. They are narrower than the fundamental (46 /an), but larger than the 7 /im predicted by lowest order 
perturbation theory. These regular profiles result from a good phase matching on axis (see Chapter III) together with 
a regular intensity dependence of the amplitude of the dipole in the plateau-cutoff transition region. Increasing the 
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intensity, this region is moved to the high density zone at the center of the jet, leading to a broadening and distortion 
of the profiles. 




FIG. 9. Radial variation of the phases of the harmonic fields corresponding to the near-field profiles shown in Fig. g: 4 xlO 14 
(dashed), 5 xlO 14 (solid) and 6 xlO 14 W/cm 2 (dot-dashed). 

In Fig. ^ we present the radial phases corresponding to these profiles. They all present a regular parabolic behavior, 
whose dependence is between 0.045r 2 and 0.048r 2 rad, r being the radial coordinate in /jm. To understand the origin 
of these curved phase fronts, let us consider the phase of the polarization at the exit of the medium. Given the low 
density, the harmonic field is obviously not mainly generated there, but this gives an estimate of what happens in the 
medium and can be directly compared to the phase of the generated harmonic field. There are two main contributions 
to the polarization phase: the first one is the Gaussian fundamental field phase multiplied by the order: 

- garc tan^)+g|(-^) , (13) 




where w{z) — woy/l + 4z 2 /b 2 and wo is the beam waist, related to the confocal parameter by b = 2itwq/\. The 
second contribution is the dipolc phase, which depends on the intensity. The harmonics are here generated in the 
plateau cutoff transition region, where the dipole phase varies linearly with intensity, with a negative slope: —77. This 
contribution can then be written as: 

-VXlfrz) = -,,_|» «p I -2( — ] + (__] , (14) 

We can here assume r << w(z), since the harmonic profiles are much narrower than the fundamental. In Chapter 3, 
we have considered the first terms of both contributions, that correspond to variations on axis. The second terms of 
these contributions, related to radial variations, both present a quadratic behavior. The term due to the focusing of 
the fundamental (Eq. (|l3|)) is in fact the phase corresponding to a Gaussian harmonic field with the same confocal 
parameter. Its radial dependence is equal to 0.032r 2 rad and is thus significantly slower than the ones observed for 
the harmonic fields. The reason for the bigger curvature of the phase fronts is thus the additional radial variation 
introduced by the dipole phase (second term in Eq. (|l4|)), which lies between 0.016r 2 and 0.023r 2 rad, depending on 
the considered intensity. 
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FIG. 10. Radial variations of the intensity and phase of the harmonic field generated at the peak intensities: 4 x 10 14 
(dashed), 5 x 10 14 (solid) and 6 x 10 14 W/cm 2 (dot-dashed), and " backpropagated" to the focus position (3mm before the jet). 



These curved phase fronts have important consequences, as we shall see below. The harmonic fields emitted at 
4 xlO 14 and 5 xlO 14 W/cm 2 present Gaussian profiles and quadratic phases at the exit of the medium, and thus 
correspond to lowest order Gaussian modes, whose characteristics (beam waist position and confocal parameter) can 
be easily calculated. We find, in both cases, a virtual focus located at the laser focus (like in the perturbative case) but 
with an extremely small size: about 1.6 /im (even smaller than the perturbative 3.7/im) The corresponding confocal 
parameters are also very small (about 1mm), so that the harmonic fields at the exit of the medium are already as if 
they were in the far-field: spherical phase front centered at the focus and linear increase of the beam size with the 
distance from the focus. Note that the radial variation of the phase of a spherical wave at a distance of z — 3.8 mm 
is: 0.046r 2 rad, close to the observed variations. 

If we calculate the angular spread of these beams far from the jet, we find similar profiles, with half angle in 1/e 2 
between 3 and 5 mrad, that can be superimposed on the "near-field" profiles if their spatial dimension is calculated 
at a distance of z — 3.8 mm. This is shown in dashed lines for the three intensities in Fig. |, where the horizontal 
scale is interpreted as the divergence (half angle at half maximum) and is simply related to the radial coordinate by: 
divergence(mrad)=0.26 xr(/im (or r — 3.8x divergence). Note that the profile corresponding to 6 xlO 14 W/cm 2 is 
more distorted by the propagation, due to its multimode structure, but its divergence is still governed by its curved 
phase front at the exit of the medium. 
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FIG. 11. Coherence degree (solid line) of the harmonic beam generated at z = 3mm and at a peak intensity of 6 xlO 14 
W/cm 2 , calculated at the exit of the medium between the central point and the outer points: for a peak pressure of (a) 3 Torr, 
(b) 150 Torr. The corresponding time integrated spatial profiles are shown in dashed lines. 

Finally, we can check the preceding conclusions concerning the virtual source size, and extend them to the case 6 
x 10 14 W/cm 2 , by "back propagating" these beams to their virtual focus. The solution of the propagation equations in 
free space for the complex conjugated field is formally equivalent to a time reversal, hence a "back propagation". Fig. 
|l0| presents the profiles and phases of the harmonic beams corresponding to the preceding intensities, calculated at 
the laser focus position {z = 0). The phases stay about constant over the extent of the profiles, indicating quasi-plane 
phase fronts. This position is thus close to the best focus for these beams. The spot sizes are very small, close to the 
above estimates, with a minimum at 6 xlO W/cm 2 with a 1.4 /an radius at 1/e 2 . 

The very small size of the virtual source of the harmonic beams is a valuable information since it gives an indication 
on the size they can be refocused to (assuming ideal refocusing without aberrations introduced by the optics, and 
no reduction factor). From this, we can estimate the achievable harmonic intensities, in view of applications. If we 
take an harmonic pulse duration of 66 fs (see Section V) and an optimized number of photons of 10 5 (obtained by 
optimizing the pressure Q), the intensity at 18.3 nm (45th harmonic) reaches about 5 xlO 8 W/cm 2 . Note that if 
argon is used as generating medium instead of neon , a m uch higher conversion efficiency is obtained for lower order 



harmonics: e.g. 10 9 photons for the 19th harmonic [ 158 1 . Since the behaviour of the dipole moment is very general, 
whatever the (sufficiently high) order or the generating gas, we can extend our results to the 19th harmonic in argon, 
and we thus find an harmonic intensity of 2 xlO 12 W/cm 2 at 43.4 nm. If a reduction factor is introduced by the 
refocusing optics, the very nice harmonic phase front should allow one to reach an even smaller focus size, and thus 
a higher intensity, provided that the surface figure of the optics is good enough. These intensities have never before 
been reached at these short wavelengths, and would open the way to nonlinear optics in the XUV region. 
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FIG. 12. Normalized spatial profile for the 45th harmonic at the maximum of the pulse for a jet position in z — —1mm 
and an intensity of 6 xlO 14 W/cm 2 . In solid line is shown the profile at the exit of the medium as a function of the radial 
coordinate, and in dashed line, the far-field profile (divergence). The fundamental near-field profile is presented in dotted line. 

So far, we have studied the beam spatial characteristics at a given time during the laser pulse. Let us now consider 
the correlation in time of the fields inside the harmonic beam, by calculating the corresponding degree of spatial 
coherence, given by Eq. (fLq). It is shown in solid line in Fig. |ll](a), calculated at the exit of the medium (z = 3.8 
mm) between the central point (on the propagation axis) and the outer points, for a peak intensity of 6 xlO 14 W/cm 2 
and a peak pressure of 3 Torr. The coherence degree is very high, and stays above 0.9 over the whole extent of the 
time integrated spatial profile, which is shown in dashed line. At this low pressure, the influence of the free electrons 
generated by the (low) ionization of the medium is negligible. The small decrease of the coherence degree can be 
understood by the slow variation over the pulse of both the phase front curvature and the spatial profile, as shown 
for three different intensities in Figs. || and|9[ At a pressure of 15 Torr, the coherence degree is a little reduced, but 
stays very high (above 0.8). However, when the pressure is increased to 150 Torr, as shown in Fig. [TT|(b), the spatial 
profile is broadened and somewhat distorted, while the coherence degree drops abruptly between 15 and 25 ^m to a 
value of 0.3. 

Indeed, the free electron dispersion is not any more negligible and results in a phase shift that is spatially dependent, 
due to the spatial distribution of the free electron density. Close to the propagation axis, the laser intensity, and thus 
the free electron density, varies slowly radially, preserving the coherence. Further away, when the laser intensity drops, 
the phase shift imparted to the harmonic field is much smaller, resulting in a decorrelation with the center of the 
beam, and therefore, a smaller coherence degree. Ionization of the generating medium is thus an important cause of 
degradation of the coherence. 

Here, we have calculated the coherence degree between the central point (on the axis) and the outer points. Note 
that, for two points taken symmetrically into the beam, the coherence degree is equal to 1, due to the revolution 
symmetry imposed to the problem. In experiments, this symmetry may be broken by different factors (gas jet and 
laser spot inhomogeneities etc.), which results in a smaller coherence degree, as recently measured by Ditmirc et at 
pif . However, our study shows that, in order to truly characterize the coherence of the harmonic beam, it is necessary 
to measure the coherence degree for non- symmetrical positions. This point makes the situation very different from 
that of a beam originating from a completely incoherent source. In the latter case, the coherence degree only depends 
on the relative position of the two points considered, i.e. their distance, not on their absolute positions into the beam, 
as it is the case here. 

In conclusion for this focus position, at low pressure the very high coherence degree at the exit of the medium 



indicates that the harmonic beam is extremely coherent, in comparison to soft X-ray lasers 153 , 32 ,146] . The coherence 
of soft X-ray lasers is often degraded by plasma fluctuations, resulting in a coherence degree equivalent to that of a 
spatially incoherent disk source with diameter of a few 100 /im [0. The associated transverse coherence length at the 
output of the X-ray laser is a few /im, whereas it reaches several tens of /jm for the harmonics. 
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C. Study of the spatial coherence: atomic jet before the focus 




FIG. 13. Radial variation of the phases of the harmonic fields generated for a jet position in z — —1mm and intensities of 
respectively: 4 xlO 14 (dashed), 5 xlO 14 (solid) and 6 xlO 14 W/cm 2 (dot-dashed). The radial variation of the polarization 
phase at the exit of the medium at an intensity of 6 xlO 14 W/cm 2 is shown in dotted line. 

The 45th harmonic spatial profile at the exit of the medium for a peak intensity of 6 xlO 14 W/cm 2 and a focusing 
1 mm after the jet is shown in solid line in Fig. [l2| It is rather distorted and exhibits an annular structure with 
three major rings. Note that very little energy is emitted on axis. As shown in Chapter III, this is the result of the 
rapid variation of the polarization phase on axis, which favors the off axis phase matching. The external radius is 27 
/im, larger than the 25 urn in 1/e 2 of the fundamental beam (dotted line). The corresponding phase is presented in 
dot-dashed line in Fig. [ll| together with phases of beams generated at 4 x 10 14 (dashed) and 5 x 10 14 W/cm 2 (solid). 
Compared to the case z = 3 mm, these phases are more irregular and vary much more quickly. Surprisingly they 
correspond to diverging phase fronts, even though the harmonics are generated by a converging beam. 

Consider the two contributions to the transverse phase of the polarization. The phase induced by the fundamental 
field, written in Eq. (|l3|), is indeed negative but varies very slowly close to the focus (quasi-plane phase front). On 
the contrary, the contribution of the atomic phase, shown in Eq. (p"4|), is associated to a very diverging phase front: 
not only the radial variation of the intensity is very rapid close to the focus, but also the harmonic is generated here 
in the plateau region where the average slope of the intensity dependence of the phase is twice as large as that in 
the cutoff region. Note that the quadratic approximation does not hold here, since we consider r ~ w{z). The total 
polarization phase at the exit of the medium and for an intensity of 6 xlO 14 W/cm 2 is shown in dotted line in Fig. 
|l3| . Its behavior is very similar to that of the phases of the harmonic beams, except for a slightly concave curvature. 

These very curved phase fronts dominate the propagation, and the calculation of the far field distribution of these 
beams gives annular profiles, as shown in dashed line in Fig. [T^, where the horizontal scale is interpreted as the 
divergence and related to the radial coordinate by: r = 1.6x divergence. The external half-angle is 15 mrad, larger 
than the 10 mrad at 1/e 2 of the fundamental. The harmonic is thus more divergent than the fundamental, which 
almost leads to a spatial separation of the two beams. 
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FIG. 14. Radial variations of the intensity (solid line) and phase (dashed line) of the harmonic field generated at z 
and 6 xlO 14 W/cm 2 and "backpropagated" to the entrance of the medium (z = —1.8mm). 
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As before, we can "back propagate" this beam and find the position of the virtual focus. The highest harmonic 
intensity is obtained in z — — 1.8 mm, which happens to be the entrance of the medium. The profile (solid line) and 
phase (dashed line) at this virtual focus are presented in Fig. 14, They look like the diffraction pattern of a ring 
with a spherical phase front. The phase shows a tt shift (i.e. a change of sign) for each profile oscillation, whereas 
the profile exhibits a narrow central peak (0.5 /im) surrounded by weak rings. A part of the energy is "lost" in these 
rings, so that the harmonic intensity reached at the central peak is barely larger than that obtained for a focusing 3 
mm before the jet. 
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FIG. 15. Coherence degree (solid line) of the harmonic beam generated at Z = —1mm and a peak intensity of 6 xlO 14 
W/cm 2 , calculated at the exit of the medium between the point located close to the center of the ring (r = 22fim) and the 
other points. The time integrated spatial profile is shown in dashed line. The peak pressure is 3 Torr. 



Let us now consider the coherence degree of this beam calculated at the exit of the medium for a peak pressure of 
3 Torr, shown in solid line in Fig. [l^. Since there is very little energy emitted on axis (see the time-integrated profile, 
shown in dashed line), we take as a reference the point located at the center of the ring (r = 22 fim). The coherence 
degree drops dramatically on both sides to about 0.3, and then oscillates around 0.2. The coherence of this beam is 
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thus very much degraded, not by the free electron dispersion (low pressure) , but by the rapid phase front fluctuations 
induced by the variation of the intensity in the laser pulse (see Fig. |l|). The dynamically induced phase can thus be 
responsible for a dramatic degradation of the spatial coherence of the harmonic beam. Note that in traditional, i.e. 
perturbative, low-order harmonic generation, where the harmonic dipole moment does not exhibit an intrinsic phase 
and depends regularly on the laser intensity, the coherence of the laser beam is simply transmitted to the harmonic 
beam. In conclusion, the dynamically induced phase plays a central role in determining the spatial characteristics 
of the generated harmonic beams, and especially their coherence properties, which are very sensitive to the focusing 
conditions. 



V. TEMPORAL AND SPECTRAL COHERENCE 
A. Influence of the jet position 
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FIG. 16. Temporal profiles of the 45th harmonic generated at 6 xlO 14 W/cm 2 , for different z positions indicated in mm in 
the captions. The laser profile is presented in dotted lines, the stars indicating the position of the inflection points. 

We have shown in the preceding chapter that the dynamically induced phase results in a dynamic phase matching 
due to the intensity distribution in the laser pulse. This will affect the temporal as well as the spectral properties of 
the generated harmonics. Let us first consider the harmonic temporal profile. Since we calculate in the slowly varying 
envelope approximation, the temporal profile is simply obtained from the response to elemental square laser pulses, 
with a Gaussian distribution of the intensities. Fig. |5j shows clearly that the temporal behavior will be very much 
dependent on the position of the jet. On both sides of the conversion efficiency curve, for example in z = — 1 mm 
and in z = 3 mm, the efficiency increases quickly with intensity. The harmonic yield will thus be high mainly for 
intensities close to the maximum of the laser temporal envelope, leading to narrow and regular harmonic temporal 
profiles. On the contrary, for intermediate z positions, phase matching can be more efficient for low intensities than 
for the maximum of the laser envelope, resulting in large and distorted harmonic profiles. 

These predictions are confirmed in Fig. [l6], which presents the 45th harmonic temporal profiles obtained at 6 x 10 14 
W/cm 2 for different z positions (ionization is not yet taken into account). When the jet is moved from z = 4 mm 
to z — 2 mm, the harmonic profile gets larger (from 36 to 125 fs FWHM) but remains regular. On the other side, 
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in z — — 1 mm, the profile is also smooth, narrow and very similar to the one obtained in z = 3 mm, with the same 
FWHM (67 fs). On the contrary, in z = 1 mm and z = 0, the profiles are quite distorted, with fast fluctuations. In 
z = 1 mm, phase matching is much more efficient at 3 xlO 14 W/cm 2 than at higher intensities (cf Fig. |5|), which 
leads to narrow peaks in the harmonic profile at half maximum of the laser envelope (shown in dotted lines). Rapid 
fluctuations also appear in the profile corresponding to the position z = 0. They are probably due to the oscillations 
of the slope in the intensity dependence of the dipole phase in the plateau region, that induce rapid changes of the 
dynamic phase matching. 
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FIG. 17. Spectral profiles of the 45th harmonic generated at 6 xlO 14 W/cm 2 , for different z positions indicated in mm in 
the captions. 

The corresponding spectral profiles are presented in Fig. [Tt]. For a focusing sufficiently before the medium, the 
spectra are quite narrow (e.g. 0.4 A FWHM at z = 3 mm). When the jet is located at z = 1 mm, the base of the 
profile gets broadened with superimposed oscillations. In z = 0, the profile is very broad with a FWHM of 4 A. 
Finally, in z = — 1 mm, the profile is again regular with a width of 2.2 A. These large variations in spectral width 
cannot be explained by variations in the intensity temporal profiles. Note for example that in z = — 1 and z = 3 mm, 
the temporal profiles are very similar with the same FWHM, whereas there is a factor of 5 between the corresponding 
spectral widths. Consequently, these broad spectra are not due to variations of the amplitude of the harmonic fields, 
but to variations of their phase. Calculations performed without taking into account their phase, i.e. assuming that 
the fields are Fourier transform limited, give indeed much narrower spectra, with widths smaller than 0.15 A. 

The origin of this phase modulation is the temporal variation of the dipole phase during the laser pulse. Let us 
model the intensity dependence of the dipole phase by a linear decrease, whose slope —r\ depends on the region of the 
spectrum: in the cutoff -q = 13.7 x 10 14 rad/(W/cm 2 ) and in the plateau r\ — 24.8 x 10 14 rad/(W/cm 2 ). The induced 
modulation of the emitted harmonic field reads: 

A$(i) = -vl(t) = -r7/ exp(-41n2(t/T) 2 ) (15) 

where r is the FWHM of the Gaussian laser envelope. This results in a modulation of the instantaneous frequency of: 

g(Ag) dl(t) 

with a corresponding broadening of the spectrum. This phenomenon is t hus very similar to self-phase modulation 



of an intense laser pulse in a Kerr medium with a negative index n-i [142]. However, the phase modulation is here 
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induced on the harmonic, whose temporal profile is different from that of the laser, and we shall see the consequences 
below. 

The rising edge of the harmonic pulse is thus shifted to the blue, and the falling edge, to the red. The extremal 
instantaneous frequencies are symmetrical on either side of the central frequency, and correspond to the inflection 
points of the laser temporal envelope: U = ±r/\/8ln2. One then finds: 

-A 2 A 2 V8hi2 , , s 

AA eKt = - — Atu ext = ±- exp -1/2 77/0 (17) 

Zttc Zttc t 

and in our conditions: AX ext ~ ±1.6 x 10 _2 77io (A). The maximal broadening of the spectrum should thus be 
observed at focus, which experiences the largest intensity variation (Iq — 6 x 10 14 W/cm 2 ) and where the intensity at 
the inflection points is still in the plateau region, ensuring the larger slope rj. This results in AA ext ~ ±2.4 A, close 
to the broadening observed for z = 0. However, this spectrum does not exhibit the characteristic grooved shape of 
the self phase modulation spectra. 



Let us recall the main features of self phase modulation [142]. Close to the inflection points, the instantaneous 
frequency varies slowly, unlike in the vicinity of the pulse maximum. The spectral density is thus maximum on the 
sides (in AA ext ) and minimal at the central frequency. On the A$(t) curve, there exist two points of the same slope, 
located on either side of the inflection point. They thus correspond to the same instantaneous frequency and interfere 
constructively or destructively depending on their relative phase. This results in a grooved spectrum, with clear 
peaks and valleys. In our case, the phase modulation is induced on the harmonic, which is generated efficiently for 



large enough intensities, located in general above the inflection points (cf Fig. 16). Consequently there will be no 
interferences in the spectra. They thus are more regular, except the one in z — 1 mm which corresponds to the only 
temporal profile sufficiently large to go past the inflection points. 
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FIG. 18. (a) Femporal dependence of the intensity (dot-dashed) and phase (solid line) of the harmonic field exiting the 
medium on axis (r = 0) for a focusing 3 mm before the jet. The variation of the polarization phase at that point is shown in 
dashed line, (b) Corresponding frequency chirps for the field (solid) and the polarization (dashed). 

In the following, we study in more details what happens in z = 3 and z = — 1 mm. In the case z = 3 mm, the 
integrated spatial profile at the exit of the medium is very regular, with a maximum on axis. The phase temporal 
variation of the harmonic field in r = has thus a large incidence on the full spectrum, and is characteristic of the 
phenomenon. Fig. |l8](a) presents the temporal variation of the intensity (dot-dashed line) and phase (solid line) of 
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the harmonic field emitted on axis at the exit of the medium. The intensity profile is quasi square, with a FWHM of 
90 fs. The phase presents a very regular behavior. Let us compare it with the polarization phase (shown in dashed 
line) calculated at this point from Eq. (|l5|). The maximum intensity there is Iq = 1.8 x 10 14 W/cm 2 (the exit of the 
medium is 3.8 mm away from the focus), which corresponds to the cutoff region (small value of 77). The two curves 
are very similar over the width of the intensity profile, confirming the origin of the phase modulation. The induced 
frequency chirps are presented in A in Fig. fl8|(b) . The spectral width corresponding to efficient harmonic emission is 
about 0.5 A, close to that of the full spectrum. Note that the chirp is almost linear during the harmonic pulse. 




FIG. 19. (a) Temporal dependence of the intensity (dot-dashed) and phase (solid line) of the harmonic field exiting the 
medium in r — 20/-im for a focusing 1 mm after the jet. The variation of the polarization phase at that point is shown in 
short-dashed line, and in a point of maximum intensity of 3.4 xlO 14 W/cm 2 , in long-dashed line, (b) Corresponding frequency 
chirps for the field (solid) and the polarization (long and short-dashed). 

In z = — 1 mm, the near-field spatial profile is annular due to an efficient phase matching off axis (cf Chapter 4). Let 
us consider the harmonic field in r = 20/jm at the exit of the medium, where the emission is efficient. The intensity 
profile, shown in dot-dashed line in Fig. [l9](a), exhibits rapid oscillations superimposed on a 90 fs wide curve. These 
fluctuations do not appear in the spatially integrated temporal profile of Fig. [16], which is very regular. This implies 
that they are compensated by opposed oscillations of the emission in neighboring points. This can be understood by 
the variation of the phase front curvature in time, which leads to fluctuations of the angle of the emission cone. 

In a given position, one then sees the radiation going back and forth, hence these oscillations. The phase of the 
harmonic field, shown in solid line, varies much more rapidly than in z = 3 mm. Since the position considered is far 
from the axis, the maximum intensity corresponds to the cutoff region : Iq = 1.7 x 10 14 W/cm 2 , and the polarization 
phase at that point (short-dashed line) varies too slowly to explain the harmonic phase behavior. A good agreement 
is obtained if we take an intensity of Io = 3.4 x 10 14 W/cm 2 , with the plateau slope 77 (long-dashed line). This implies 
that the harmonic radiation which exits the medium at that point has mainly been generated in a deeper region, close 
to the axis. The phase front curvature makes it diverge and exit at that position. The frequency chirp, presented in 
Fig. |l^(b), shows the phase sudden changes of slope, reflecting the intensity fluctuations. The corresponding spectral 
width is about 2.2 A, close to that of the full spectrum. 

The origin of the difference by a factor 5 of the spectral widths in z — 3 and z = — 1 mm appears thus clearly: 
in z = 3 mm, the harmonic emission corresponds to the cutoff region, hence a small slope 77 and a low intensity. In 
z = — 1 mm, it is mainly in the plateau, with a slope twice as large, and for intensities at least two times larger. 
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The spatial and spectral coherence properties are thus closely linked, and are both governed by the variations of the 
atomic phase. 

B. Influence of the ionization 
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FIG. 20. Temporal (a) and spectral (b) profiles for a focusing 1mm after the jet at 6 xlO 14 W/cm 2 : without ionization 
(solid), with ionization at a pressure of 3 Torr (short-dashed) and 15 Torr (long-dashed) of neon. In dot-dashed line is shown 
the profile corresponding to 15 Torr when the defocusing of the laser is neglected. 

When the jet is placed several mm away from the focus, the intensity in the medium is sufficiently low to neglect 
the effects of ionization. In z = 3 mm for example, the profiles are only slightly modified. We will thus concentrate 
on the case z = —1 mm, and try to separate the different effects: depletion of the atomic population, defocusing of 
the laser and phase mismatch both due to the generated free electrons. 

Fig. |2^(a) presents the distortions of the temporal profile induced by the different effects. In short-dashed line is 
shown the profile corresponding to a peak pressure of 3 Torr. This low pressure ensures that the free electron density 
is low enough to neglect its effects. This profile thus characterizes the incidence of the depletion of the medium on 
harmonic generation (the susceptibility of the generated ions is supposed to be negligible compared to the atoms) . At 
the center of the jet, ionization reaches 67% at the end of the laser pulse. The resulting profile is clearly asymmetric, 
with a rising edge hardly modified and a decrease of the efficiency on the falling edge. The FWHM is slightly decreased, 
from 67 to 63 fs. 

At a pressure of 15 Torr, the effects of the free electrons are not any more negligible. We have first made the 
calculations by neglecting the defocusing of the laser: we suppose a Gaussian propagation and we only solve the 
propagation equation for the harmonic field in the ionized medium. The temporal profile (dot-dashed) is even more 
distorted on the falling edge: the width decreases to 58 fs. This is due to the phase mismatch induced by the free 
electron dispersion, that reaches Ak ~ 22 mm -1 at the center of the jet and at the maximum of the pulse. 

The study of the propagation of the laser in the ionizing medium shows clearly an effect of defocusing: at the exit 
of the medium (0.2mm from the focus), the maximum intensity is reduced by 17%, to 5 xlO 14 W/cm 2 . The reduced 
intensity experienced by the medium results in a smaller ionization degree: 54% at the center of the jet and at the 
end of the pulse (vs 67%). A consequence of the defocusing is thus to reduce the amplitude of the preceding effects: 
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less depletion and smaller free electron density. However, the corresponding temporal profile (long-dashed) is the 
most distorted one. The efficiency drops abruptly after the maximum of the pulse, leading to a width of 46 fs. Two 
elements can explain this decrease in efficiency. First, the reduced intensity induces a decrease of the polarization 
amplitude. Second, the defocusing changes the intensity distribution in the medium compared to the Gaussian. The 
atomic phase distribution is thus modified, which may result in a less efficient phase matching. 

In Fig. ^o|(b) we show the incidence of these phenomena on the spectrum. A striking feature is that, whatever the 
effect considered, the distortions of the spectrum simply follow the distortions of the temporal profile. The blue side 
of the spectrum is weakly affected, whereas its red side is progressively cut off. The FWHM decreases first from 2.2 to 
1.9 A due to the depletion (3 Torr), then to 1.6 A due to the phase mismatch induced by the free electrons at 15 Torr, 
and finally to 1.5 A when the defocusing of the laser is taken into account. This is a clear illustration of the phase 
modulation phenomenon. Since the rising edge of the harmonic pulse is weakly affected by ionization, it is also the 
case for the blue side of the spectrum it is associated to. On the contrary the falling edge, and consequently the red 
side of the spectrum, are strongly distorted. This results in a blueshift of the central frequency of the spectrum. Note 
that this phenomenon is differe nt from t he blueshift induced by the temporal variation of the free electron density, 
and thus of the refractive index |15£,12£|. In our conditions, the latter effect shifts the fundamental spectrum by 2.4 
A, which results for the 45th harmonic in a shift of only 0.05 A. 



C. Consequences of the phase modulation 
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FIG. 21. Compressed temporal profiles obtained at 6 xlO 14 W/cm 2 and at a pressure of 3 Torr (solid) and 15 Torr 
(dot-dashed) of neon: (a) in z = 3mm, (b) in z = —1mm. In dashed lines are shown the corresponding emission profiles 
before compression when ionization is neglected. 

The phase modulation studied below can be used to control the temporal and spectral properties of the harmonic 
radiation. In the time domain, the regular close-to-quadratic phase, and the corresponding linear chirp, lead to the 
possibility of compressing the harmonic pulse with a pair of gratings, as is done in CPA lasers. Schafer and Kulander 
studied the single atom response to very short laser pulses (27 fs), and showed that it should be possible to compress 



the harmonic pulse to durations below the fundamental period (2.7 fs) [ 139 1 . We will discuss the short duration 



problem in the next section, and we concentrate here on the study of the compression of the macroscopic harmonic 
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beam. 

To simulate the compression, we subtracted a mean quadratic phase from the phase of the spectrum emitted at each 
point at the exit of the medium. The inverse Fourier transform followed by the spatial averaging give the compressed 
macroscopic temporal profile. The result of the compression in z = 3 mm for a peak pressure of 3 Torr is presented 
in solid line in Fig. |l](a). The profile is very narrow, with a FWHM of 10 fs (vs. 66 fs before compression, dashed 
line). This is close to the Fourier transform limit. The short duration is maintained when the pressure is increased to 
15 Torr (dot-dashed line). In z = —1 mm, the compressed profile for a pressure of 3 Torr is even narrower, as shown 
in solid line in Fig. pT](b) : 6 fs vs. 67 fs before compression. 

Note that on the basis of the spectral width which is more than three times larger than in z — 3 mm, we would 
have expected a smaller duration. However, this optimal compression is not reached due to temporal fluctuations 
in the phase (see Fig. |l9|) and to spatial inhomogeneities of the phase modulation. Phase matching has thus to be 
considered carefully in order to get the shortest harmonic pulses. At a pressure of 15 Torr, not only is the spectral 
width reduced (due to defocusing of the laser, see preceding section), but also the free electron dispersion introduces a 
temporally-varying phase shift on the harmonic emission. This results in a broader temporal profile of 13 fs FWHM, 
shown in dot-dashed line. In conclusion, we have shown that it is possible to compress the macroscopic harmonic 
beam generated by a 150 fs laser to durations below 10 fs. 




FIG. 22. Spectral profiles of the 45th harmonic generated with a fundamental without chirp (solid), chirped positively 
(dashed) or negatively (dot-dashed): (a) in z — 3mm, (b) in z — —1mm. Ionization is not taken into account. 

In the frequency domain, it should be possible to compensate for (or increase) the dynamically induced harmonic 
chirp by introducing an appropriately designed frequency chirp on the fundamental beam. We have calculated the 
harmonic spectrum generated by a laser pulse of same duration (150 fs) but presenting a quadratic phase in time, 
such that the fundamental spectrum is broadened to 32 nm, corresponding to a Fourier transform limit of 25 fs. 
Fig. p2|(a) presents the spectra obtained for positive and negative chirps in the case of a focusing 3 mm before the 
jet. In the preceding section, we have seen that for this position, the dynamically induced harmonic chirp, which 
is negative, is limited (cf Fig. |lj), leading to a quite narrow spectrum (0.4 A). It is much smaller than the chirp 
induced by the fundamental on the harmonic (recall that the fundamental phase is multiplied by the nonlinear order 
in the polarization). In these conditions, we thus get a broadening of the harmonic spectrum whatever the sign of the 
fundamental chirp. 

However, when the fundamental chirp is negative (blue before the red), it adds to the negative dynamic chirp, 
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leading to a 3.4 A width (dot-dashed line). On the other hand, the positive fundamental chirp gets subtracted from 
the dynamic chirp, giving a 2.5 A width (dashed line). Note that the fundamental chirp could be adjusted to reduce 
the width of the spectrum for this particular position. In z = — 1 mm, the dynamic harmonic chirp is larger resulting in 
a 2.2 A wide spectrum (solid line in Fig. ^2|(b)). It is in fact of the same order as the one induced by the fundamental 
chirp. When both add (negative fundamental chirp), this results in an extremely wide spectrum (5.6 A, dot-dashed 
line). Conversely, when they get subtracted (positive fundamental chirp), they compensate each other, leading to a 
very narrow spectrum of 0.4 A width (dashed line). 




FIG. 23. Spectral profiles of the 45th harmonic generated in 15 Torr of neon at 6 xlO 14 W/cm 2 , in z = —1mm and with 
a fundamental chirped positively with (resp. without) ionization: short-dashed (solid), or negatively with (resp. without) 
ionization: dotted (dot-dashed). 

So far, we have not yet taken ionization into account, which may modify the results in z = — 1 mm at a sufficiently 
high pressure. At a pressure of 3 Torr, only depletion plays a role, and its influence is quite limited (see Fig. pp| ). 
The consequences of ionization for a pressure of 15 Torr are shown in Fig. ^3| The wide spectrum (negative chirp) 
is reduced to 4.0 A (dotted line) by the same effect that was described in the preceding section: loss of the red part 
of the spectrum due to a reduced efficiency on the falling edge of the pulse. On the contrary, the "compensated" 
spectrum (positive chirp) is broadened to 1.1 A (short-dashed). This is certainly caused by phase fluctuations induced 
by the free electron dispersion. 

In conclusion, the regular behavior of the dynamically induced harmonic chirp can be used to control the temporal as 
well as the spectral properties of the harmonic emission. Concerning the latter, we find a qualitative agreement with the 
experimental results of Zhou et al. 168 1: for negative fundamental chirps, they observe a very pronounced broadening 
whereas for positive chirps, the peaks remain relatively narrow. However, the experiments were performed in argon at 
a relatively high intensity, thus with a higher degree of ionization than considered here. Harmonic generation occurs 
then only on the rising edge of the laser pulse, resulting in significant redshifting or blueshifting (depending on the 



D. Influence of nonadiabatic phenomena 



So far, we have considered harmonic generation by laser pulses of 150 fs duration. The recent advent of lasers of 
duration 25 fs and less, sufficiently intense to generate harmonics, motivates the investigation of shorter pulse effects. 
Except for the influence of ionization which depends directly on time, our results are simply scalable to shorter 
durations as long as the approximations made still apply: for example, going from 150 to 25 fs laser pulses (while 
keeping the same peak intensity) would result in a shortening of a factor 6 of the harmonic temporal profile and a 
corresponding broadening of the same factor of the spectrum. What about our approximations? 

The slowly varying envelope approximation for the propagation of the harmonic field is valid for very short laser 
pulses since the harmonic period is much less than 1 fs. However, the very broad spectra of the generated harmonics 
may overlap and interfere, making the propagation calculations more difficult. But the main problem is to take 
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into account the nonadiabatic response of the nonlinear polarization to the fast driving field. In our calculations, we 
consider that the harmonic dipole moment follows adiabatically (or reacts instantaneously to) the intensity distribution 
of the laser pulse. In particular, we consider that the laser field amplitude does not vary significantly over the optical 
period. For the very short pulses mentionned above, this may not be true and result in a distorted atomic response. 

Such a nonadiabatic effect wa s inv oked by Schafer and Kulander to explain the phase characteristics of harmonics 
generated by a 27 fs laser pulse [ 139 1 . They calculated the single argon atom response by integrating numerically the 
time dependent Schrodinger equation in the single active electron approximation. The plateau region of the spectrum 
is highly structured, the expected odd harmonics being hard to distinguish. However, harmonics at the end of the 
plateau are found very broad but distinct, exhibiting a nice quadratic phase. They interprete this pheno men on in 
terms of the change in laser intensity during the laser period that alters the harmonic generation process [162]. An 
electron that enters the continuum while the laser intensity is increasing experiences an additional acceleration before 
returning. This produces a blucshift on the rising edge of the laser envelope. Conversely, electrons ionized after the 
peak of the pulse are decelerated, returning later, leading to a redshift of the spectrum. 




FIG. 24. Temporal envelope (lower curves) and phase (upper curves) of the 49th harmonic generated in argon by a 27 fs 
laser pulse at 3 xlO 14 W/cm 2 : adiabatic (short-dashed) and nonadiabatic (solid) calculations. The temporal profiles after 
compression are shown in dot-dashed (adiabatic) and long-dashed (nonadiabatic). 
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FIG. 25. Spectral envelope (lower curves) and phase (upper curves) of the 49th harmonic generated in argon by a 27 fs 
laser pulse at 3 xlO 14 W/cm 2 : adiabatic (short-dashed) and nonadiabatic (solid) calculations. The adiabatic phase minus its 
quadratic component is shown in dot-dashed line. 

In order to investigate these effects, we calculated the nonadiabatic single atom response to the fast driving field 
in the strong field approximation and compared it with the adiabatic result. We use the same conditions as in Ref. 



[ 139 1 , i.e. an argon atom interacting with a 27 fs, 810 nm pulse with a peak intensity of 3 xlO 14 W/cm 2 . In these 
calculations, we neglect ionization since the considered intensity is (slightly) below saturation for this small pulse 
duration. 

The 49th harmonic temporal envelopes and phases are presented in Fig. |24|. The adiabatic profile (short-dashed) is 
very regular with a 7.6 fs FWHM, while the corresponding phase (upper short-dashed) exhibits a close-to-quadratic 
temporal dependence. Indeed, at the considered intensity, the 49th harmonic is in the cutoff region, where the dipolc 
varies regularly with intensity. In particular, the phase depends linearly with intensity, which leads to a temporal 
quadratic phase close to the top of the Gaussian laser pulse. The nonadiabatic profile (solid) exhibits the same width 
but is delayed compared to the adiabatic by roughly 1.3 fs, whereas the phase is close to the adiabatic calculation. 
This shift in time can be related to the physics of the process. In the two-step model, there is one main electron 
trajectory for the harmonic emission in the cutoff region. The corresponding return time, i.e. the time between 
tunneling and recombination, is about half a period. Since the amplitude of the emission is mainly determined by the 
ionization probability, this explains the delayed response of the atomic dipole. 

The spectral envelopes and phases are presented in Fig. ^5| The adiabatic profile is very broad (1.3 eV FWHM), 
even more than the nonadiabatic (1.2 eV). The latter is shifted to the red, which is a consequence of the shift in 
time of the corresponding temporal profile. Indeed, the blue spectral components are associated to the decreasing 
part of the temporal quadratic phase, and are thus affected by the delay of the rising edge of the temporal envelope. 
Furthermore this positive temporal quadratic phase results in a negative spectral quadratic phase, as shown in Fig. 
p5| (upper lines). The adiabatic phase minus its quadratic component, shown in dot-dashed line, is almost constant. 
This operation can be performed experimentally with a pair of gratings, leading to the compression of the harmonic 
temporal profile. They are presented in Fig. |i| The adiabatic (dot-dashed) and nonadiabatic (long-dashed) profiles 
are similar with a 1.7 fs FWHM. The compression by more than a factor 4 res ults in a corresponding increase in peak 
intensity. All these profiles are very similar to the ones reported in Ref. |139[ |. 

In conclusion, the main features in the cutoff region such as broad spectrum and quadratic phase are very well 
reproduced by an adiabatic calculation. They are thus the consequence of the intensity dependence of the phase of 
the adiabatic dipole, and not of a nonadiabatic effect due to the change in intensity during the laser period. In fact, 
harmonics in the cutoff are created close to the peak of the laser pulse, where the intensity varies slowly. Nonadiabatic 
phenomena may play a role in the plateau region, where the harmonics are also generated on the edges of the laser 
envelope. In the cutoff, the main nonadiabatic effect seems thus to be the shift in time of the atomic response. Note 
that if we take a sine laser field instead of a cosine, the harmonic temporal profile is slightly modified, exhibiting 
a double peak structure for the same width. This is an indication that we are close to the limit of validity of the 
adiabatic approximation. 

Another important element is the width of the spectrum. It spans one laser photon energy (1.5 eV) on either side of 
the harmonic frequency. This means that for shorter pulse durations, the spectra of adjacent harmonics will overlap 
and interfere, leading to a perturbed behaviour also in the cutoff region. Note that a nonadiabatic effect has been 
invoked by Christov et al. to explain the delayed ionization with short laser pulses (25 fs) and the corresponding 
increase in harmonic plateau extension |3q] . 



VI. FUTURE APPLICATIONS 



High order harmonics are currently used in a number of applications such as atomic and molecular spectroscopy, 
solid state physics... (for a review, see Ref. f9§|l). We focus here on potential applications involving their coherence 
properties. Indeed, the main message of the studies presented in Chapters 3-5 is that harmonics can be generated in 
the form of short, but spatially and temporally coherent pulses. Moreover, the properties of such pulses can be to 
a great extent and on many aspects controlled. Therefore, one can foresee three very interesting classes of possible 
applications of harmonic pulses: applications in nonlinear optics in the XUV, that employ the high refocused harmonic 
intensities, applications in interferometry, that employ spatial coherence properties, and applications in attosecond 
physics, that employ the phase modulation and other phase properties to reduce the duration of harmonic pulses. In 
this Chapter we present a somewhat speculative discussion of these three classes of applications, since none of them 
has been realized so far. 
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A. Nonlinear optics in the XUV 



In Chapter 4, we have shown that the good wavefront quality of the harmonic beam could allow us to refocus it 
quite easily to very small spots of a few /im in diameter using simply a spherical mirror. Recently, the Saclay group has 
confirmed this prediction by measuring the size of the focal spot, and its dependence on the medium characteristics 
[§8j. Even better focusing could be achieved by using a Schwarzschild objective [Q. 

If the short pulse character of the emission is preserved, which implies the use of multilayer mirrors to select and 
refocus the harmonic beam, unprecedented high intensities could be reached in the XUV region. For instance, if argon 
is used as gene rating medium, an optimized number of 10 10 photons could be generated around 30 eV. Schafer and 
Kulandcr [139| found that the conversion efficiency at the single atom level increases roughly as quickly as the laser 
pulse duration decreases, resulting in the emission of about the same number of harmonic photons. Using the intense 
25 fs infrared lasers currently available, an harmonic pulse duration of less than 10 fs could be obtained (note that 
the generation of sub- femtosecond pulses is under consideration, see below). The refocused harmonic intensity would 
thus reach 10 14 W/cm 2 , enough to induce multiphoton processes. 

Some attempts have already been done to observe non-resonant two-photon ionization of rare gases 
They are for the moment unsuccessful!, except for moderate-order harmonics ]73| , but the rapid progress in short 
pulse laser technology (repetition rate, pulse duration and energy) will improve the performances of the harmonic 
sources and should make possible the observation of such processes in the near future. 



B. Interferometry with harmonics 



Interferometry with harmonics is particularly interesting because of their high frequency and short duration. Co- 
herent harmonics of high frequencies would pass many interesting media without significant absorbtion. Of particular 
interest are the dense laser-induced plasmas. Short wavelengths are much less sensitive than visible ones to refraction 
by the large density gradients involved, and correspond moreover to higher critical densities. Furthermore, due to their 
short duration, the harmonics may probe dynamics of many systems in a quasi-instantaneous way. Interferometry with 
harmonics thus appear as a powerful diagnosis for different medias, and in particular plasmas. The broad harmonic 
bandwith (correlated to the short pulse duration) makes rather difficult to use an amplitude division interferometer 
(Michelson's type) at the harmonic frequency, but the good spatial coherence is well suited for wavefront division 
interferometers (Young's type). Two kinds of possible interfero metr y experiments are discussed by experimentalists: 



in one kind, two correlated harmonic beams are used (like in 1 167 ), out of which one may, and the other does not 



pass through the medium to be diagnosed; both beams interfere thereafter. Realization of such experiment may pose, 
however, difficulties in achieving sufficient spatial and temporal overlap of the two harmonic pulses. 

The other kind of experiment involves only one harmonic beam which would pass through a region of spatially 
non-uniform index of refraction, and produce a self-interference pattern thereafter. Numerous applications of this last 
scheme in solid state and plasma diagnostics are currently being considered by experimentalists in Saclay. 

Note finally that the harmonic spatial coherence could be used in holographic or phase-contrast imaging. In 
particular, the harmonic generation spectrum has now reached the water window (between the K-edges of carbon at 



4.4 nm and oxygen at 2.3 nm [144 ,B3). This region is of particular interest for biological applications since it provides 



the best contrast for hydrated carbonated structures such as living cells. 



C. Attosecond physics 



The shortest pulses achieved today (5 fs in the infrared jX8|) are limited by the long period of the radiation. Shorter 
wavelengths are required for further pulse shortening. High-order harmonics are the most promising way to generate 
sub femtosecond, i.e. attosecond (as), pulses. Different schemes have been proposed to reach these extremely short 
durations. The first one deals with the relative phase of the harmonics generated in the plateau region |3^|56j. If 
they were emitted in phase, the corresponding temporal profile would consist of a train of pulses separated by half the 
laser period and of duration in the attosecond range. There is a clear analogy here with mode-locked lasers. However, 
early calculations of the single atom response showed that the high harmonics were, in general, not in phase, due to 
the interference of various energetically allowed electronic trajectories leading to the harmonic emission. Recently, 
Antoine et al. revived this proposal by showing that the propagation in the atomic medium could select one of these 
trajectories, resulting in the macroscopic emission of a train of » 200 attosecond pulses 0. 

In principle, the proposals based on phase locking require filtering of the phase locked harmonics from the entire 
spectrum; in particular they require suppression of the low harmonics from the fall-off region of the spectrum. It is 
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worth mentioning that trains of attosecond pulses could also be generated in laser-induced plasmas from the surface 
of a solid target. H armonic generation from solid targets has been intensively studied in the recent years both 
experimentally ]3l| , |74 ,101,lll|, and theo retically [33 125 IOC], but only very recently have the temporal aspects of 
such a generation process been studied [130,131]. Using a simple model of the oscillating plasma surface, as well 
as a ID "Particle-In-Cell" code, these authors have shown that the reflected signal obtained in the case of normal 
incidence of the driving laser has the form of a train of ultrashort pulses of duration ~ 20 times shorter than the 
optical period (i.e. ~ 100 attoseconds for the case of Ti- Sapphire laser). This is quite surprising because the harmonic 
spectrum is monotonically decreasing and does not exhibit any plateau, and moreover, the production of those pulses 
does not require any filtering. 

Furthermore, Corkum et al. |59 67 proposed a way of generating a single pulse by using the high sensitivity of 
harmonic generation in gases to the laser polarization j2^^]. By creating a laser pulse whose polarization is linear 
only during a short time, close to a laser period, the emission could be limited to this interval. Propagation calculations 
validated this idea p^ j. In their original proposal, Corkum et al. suggested to use two cross polarized fundamental 
beams of slightly different frequencies in order to obtain the desired form o f the time dependent polarization. Very 
recently, the Lund group has realized this idea in a somewhat different way [16C], by using a flat birefringent crystal 
and a short input pulse that is chirped in frequency and polarized at 45 degrees relative to the optical axis of the 
crystal. The input chirped pulse in such a system is split into two orthogonally polarized chirped pulses delayed in 
time. When the phase difference between the two polarizations is zero, linear polarization is produced. This enables 
a single ultrashort time window to be selected for efficient harmonic generation. Moreover, two or more ultrashort 
time windows with controllable delay may be selected. 

Another scheme of ultrashort pulse generation has been described in the preceding section, and concerns a single 
harmonic. The compression of the dynamically induced frequency chirp of an harmonic located at the end of the 
plateau and generated by a 27 fs laser pulse can result in a pulse of 1.7 fs duration. Moreover, Schafer and Kulander 
proposed to compress three adjacent harmonics to generate sa 400 as pulses [ 139 1 . This is made possible by the fact 
that harmonics located close to the cutoff exhibit approximately the same frequency chirp. 

Finally, the last scheme uses very short duration laser pulses to shorten even more the harmonic emission, thanks 
to the rapid ionization of the medium (less than one cycle). As already mentioned, with 5 fs infrared pulses it has 
beco me recently possible to generate soft XUV radiation extending to the K-edge of carbon at 4.4 nm (water window) 
[ 144 1 . Simulations show that the generated pulses could be of sub-500 attosecond duration J37j ]. 

Although the problem of generating attosecond pulses is already difficult, it is even more challenging to consider their 
detection, and applications. So far, the literature on those subjects hardly exist, but theoretical studies have begun. 
Corkum and his collaborators |h| discuss possible applications of subfemtosecond pulses in physical chemistry. These 
authors discuss, for instance, Coulomb explosion measurements for molecules in strong laser fields; such measurements 
allow to measure time-dependent molecular structures. The use of attosecond pulses for such measurments would 
allow to "freeze" the nuclear motion for even highly charged molecular ions by their inertia for the duration of the 
pulse, and thus achieve much better temporal resolution than it is possible so far. 

An experimental detection of the train of attosecond pulses is equivalent to the detection of phase locking. The 
latter can be studied for instance in the process of two-photon ionization of Helium atoms by a filtered harmonic 
signal consisting of 5-10 harmonics. The ionization yield in such a process should depend on the relative phases 
between the harmonics 0]. Maquet and his collaborators in a series of papers has investigated the possibility of 
determining the relative phases between the harmonics and the fundamental field from the above threshold ionization 



(ATI) spectra. In Refs. [156,147] ATI in a two color field (fundamental plus one high harmonic frequency) was 



cansidered. The spectra in that case may be strongly affected by quantum interferences between different ionization 
and electron redistribution channels, and depend on the relative phase between the two fields. In further publication 
[ 157 1 the authors considered ATI in a mixed multicolor field consiting of the fundamental and several high harmonic 
components. The detailed study of photoionisation spectra provides in this case more, or less direct information on 
the phase differences between successive harmonics. Comparison with experimental results M suggests that indeed 
harmonics generated in rare gases are phase locked. 

Another possibility of detection and application of attosecond pulses would be to develop the time resolved attosec- 
ond spectroscopy (TRAS). This could be done in analogy with femtosecond spectroscopy, using a pump and probe 
attosecond pulses, for example to excite and then ionize a coherent sub-Rydberg wave packet in an atom. The result 
of such pump-probe process should depend on the time delay between the pump and the probe. 

In the review [jl0| we suggested that TRAS could also be realized by mixing an attosecond pulse train (the probe) 
with the strong fundamental laser beam (the pump). In particular, we have considered harmonic generation by 
such multicolored fields, expecting that the resulting harmonic signal would depend on the time delay (modulo 
one fundamental period) between the train and the pump. Unfortunately, as we already pointed out in jlQ] , the 
estimates presented in that paper, based on SFA, are too optimistic. We have checked, using TDSE method, that the 
dependence of the harmonic signal on the time delay in such process is much weaker in the experimentally accessible 
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regime of parameters. Nevertheless, the scheme allows at least to study temporal interference between the train of 
pulses focused in the medium (probe), and the train of pulses generated in the medium by the fundamental laser 
(pump). This interference can be particularly efficient due to phase matching effects, such as trajectory selection via 
propagation JtJ. 

Other processes than harmonic generation, can be considered in order to make TRAS with attosecond pulse trains, 
or separated attosecond pulses possible. In Rcf. [148 ] the authors consider ATI process by a fundamental laser pulse 
that leads to the resonant excitation of dynamically Stark-shifted Rydberg states. Short harmonic pulse (of duration 
much shorter than the fundamental one) is applied during that process, and its delay with respect to the front of the 
fundamental pulse is controlled. As a result, the highly structured ATI spectra depend on the delay of the probe, and 
provide thus the first example of time-resolved spectrocopy of dynamically induced resonances. The quest for other 
forms of TRAS remains, in our opinion, one of the most interesting challenges of super-intense laser atom physics. 



VII. CONCLUSION 



We have used a theory of high order harmonic generation by low frequency laser fields in the Strong Field Approxi- 
mation to study the spatial and temporal coherence properties of the harmonic radiation. We show that the intensity 
dependence of the atomic dipole phase predicted by this model plays a central role in the way phase matching is 
achieved in the medium. In particular, it introduces an asymmetry relative to the position of the jet compared to the 
laser focus. When the laser is focused sufficiently before the jet, phase matching is optimized on axis, leading to very 
regular Gaussian emission profiles. When it is focused in or after the medium, harmonic generation is prevented on 
axis but can be efficient off axis resulting in annular profiles. Moreover, the dipole phase induces a big curvature on 
the generated harmonic phase front, leading to quite diverging harmonic beams. Consequently the virtual source of 
the emission presents a very small size. Temporal fluctuations of the harmonic phase front can degrade the spatial 
coherence of the radiation, in particular for a focusing after the jet. However, before the jet, the spatial coherence 
can be very high indicating that high order harmonics could be a useful coherent source in the XUV. 

Induced by the intensity distribution in the laser pulse, the dynamic phase matching results in very regular, narrow 
temporal profiles for focus positions sufficiently before or after the jet. At intermediate positions, the harmonic 
profiles are distorted and almost as large as the fundamental. The corresponding spectra are very broad and far from 
the Fourier transform limit. This is due to the intensity dependence of the dipole phase, which results in a phase 
modulation of the emitted harmonic field. This phase modulation is all the more important as the peak intensity is 
high and the harmonic emission is in the plateau region, where the slope of the intensity dependence is the larger. 
The harmonic beam thus presents a negative chirp, with a rising edge shifted to the blue and a falling edge, to the 
red. In our conditions, ionization is quite limited and affects mainly the falling edge, and thus the red part of the 
spectrum. 

The regular phase modulation can be used to control the temporal and spectral properties of the harmonic radiation. 
Temporally, the chirped pulse could be recompressed by a pair of gratings to very small durations. Spectrally, an 
appropriately designed chirp on the fundamental can compensate (or increase, depending on its sign) the harmonic 
dynamic chirp, resulting in a narrow (or extremely broad) spectrum. Finally, we show that cutoff harmonics generated 
by pulses as short as 27 fs can be described by an adiabatic calculation, the nonadiabatic phenomena being negligible. 

In conclusion, the spatial and temporal coherence properties of high order harmonics are closely related, and depend 
strongly on the atomic dipole phase. However, the regular behavior of the latter makes it possible to control them by 
using relatively simple means. 

These unprecedented coherence properties in the XUV region open the way to a new class of potential applications 
which are attracting more and more interest. The good wavefront quality, that allows focusing to very small focal 
spots, together with the short pulse duration could result in the high focused intensities necessary to perform nonlinear 
optics in the XUV. The high spatial coherence makes it possible to realize XUV interferometry experiments, such 
as ultrashort plasmas diagnostics. Finally, the harmonic phase characteristics in time could be used to generate 
attosecond pulses, the shortest ever, opening the way to "attosecond physics". 

We acknowledge fruitful discussions with P. Agostini, Th. Auguste, K. Burnett, B. Carre, P. B. Corkum, T. Ditmire, 
M. Dorr, M. Gaarde, M. Y. Ivanov, Ch. Joachain, H. C. Kapteyn, K. C. Kulander, A. Maquet, P. Monot, M. M. 
Murname, M. D. Perry, A. Sanpera, K. J. Schafer, R. Ta'ieb, V. Veniard, C.-G. Wahlstrom, and J. Watson. 
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